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We study the general structure of models for structure formation, with applications to the reverse 
engineering of the model from observations. Through a careful accounting of the degrees of freedom 
in covariant gravitational instability theory, we show that the evolution of structure is completely 
specified by the stress history of the dark sector. The study of smooth, entropic, sonic, scalar 
anisotropic, vector anisotropic, and tensor anisotropic stresses reveals the origin, robustness, and 
uniqueness of specific model phenomenology. We construct useful and illustrative analytic solutions 
that cover cases with multiple species of differing equations of state relevant to the current generation 
of models, especially those with effectively smooth components. We present a simple case study of 
models with phenomenologies similar to that of a ACDM model to highlight reverse-engineering 
issues. A critical-density universe dominated by a single type of dark matter with the appropriate 
stress history can mimic a ACDM model exactly. 
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I. INTRODUCTION 

How does one reverse engineer a model for structure 
formation from observed phenomena? How unique is 
such an inversion? How robust are the phenomenolog- 
ical distinctions between broad classes of models? With 
the wealth of high-precision cosmological data expected 
in the near future from the cosmic microwave background 
(CMB), galaxy surveys, and the high redshift universe, 
the simple ab initio models for structure formation cur- 
rently considered may be ruled out, forcing us to con- 
front these difficult issues. In this paper, we take the 
first steps toward answering this question by examining 
from a general standpoint what makes a model for struc- 
ture formation behave as it does in linear perturbation 
theory. 

A model for structure formation is completely speci- 
fied by its initial conditions and the full temporal and 
spatial behavior of the stresses in its dark sector. The 
dark sector contains the elements in the model that do 
not interact with the photons at any observable redshift. 
It can include, but is not limited to, cold dark matter 
(CDM), neutrinos, and cosmological defects. 

Unfortunately, the stress history of the dark sector is 
by definition not directly observable. Its effects come fil- 
tered through gravity as mediated by metric fluctuations. 
The translation of metric fluctuations into observables in 
the CMB and evolution of structure is well understood. 
Therefore, the main hurdle in the task of reconstructing a 
model from observations is to understand how stress his- 
tories translate into metric fluctuations and vice versa. 

Our general philosophy here is to start from elements 
of the cosmological model that will likely survive the on- 
slaught of data: general relativity and a universe whose 
deviations from homogeneity and isotropy are initially 
small. We proceed down the theory pipeline to the ex- 
isting models of structure formation, making explicit the 
places where assumptions are made and hence could be 
altered. Where possible, we provide analytic solutions 



and approximations that highlight certain generic be- 
havior and phenomena. These solutions are useful for 
describing the behavior of the existing models and are in 
most cases new or substantially more general than those 
found in the literature. In particular, we derive master 
solutions for models which contain multiple components 
with arbitrary equations of state and smooth, entropic, 
anisotropic, and sonic stresses. 

We begin in §|| with an overview of the basic elements 
of a structure formation theory and their traditional clas- 
sification in terms of their initial conditions, perturbation 
type, and clustering properties. We then present a con- 
cise but general treatment of linear perturbation theory 

Although these are well- 
i, our treatment has sev- 



IV 



in § [II and gauge issues in 
studied subjects (see e.g. j 
eral pedagogical and practical virtues. It keeps careful 
track of the degrees of freedom available to structure for- 
mation models and hence provides a unified treatment 
applicable to all models, including those containing ex- 
otic matter like scalar fields or cosmological defects. We 
also explicitly maintain general covariance such that the 
equations apply, and may be easily specialized, to any 
choice of coordinates or gauge. In §|v|, we define general 
classes of stress perturbations and present an overview of 
their conversion into observables. 

The remainder of the paper deals with stress histo- 
ries on a case by case basis. The simplest case involves 
smooth stresses, and we present detailed analytic solu- 
tions in § VI that apply to a wide range of models — from 



simple cosmological constant models (ACDM) to massive 
neutrino and scalar field models. Pure aniso tropic, en- 
tropic, an d son ic stresses are treated in § |VIl| and mixed 
cases in § VIII . To highlight reverse-engineering issues, 
we study single-component, critical-density dark mat- 
ter models with phenomenologies that mimic the ACDM 
model in §IX. We conclude in §[x] by re-examining the 
traditional classification scheme of §|l| in light of the phe- 
nomenological distinctions uncovered in this work. 
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II. CLASSIFICATION OF THEORIES 

A. Initial Conditions 

Perhaps the most fundamental difference between 
models for structure formation lies with their initial con- 
ditions. Currently, inflation is the only known means of 
laying down large-scale density or curvature perturba- 
tions in the early universe. Indeed, inflation in the more 
general sense of a period of superluminal expansion is 
required for the causal generation of large-scale power. 
It provides a means of connecting parts of the universe 
that are currently space-like separated, i.e. outside the 
current particle horizon. Models with initial curvature 
perturbations are usually called "adiabatic" models. 

All other causal models begin with no density or cur- 
vature fluctuations on large scales and are hence called 
"isocurvature" models. In these models, stress gradients 
causally move matter around inside the horizon to form 
large-scale structure. 

The generation mechanism is also responsible for de- 
termining the spectrum and statistics of the fluctuations. 
The simplest inflationary models predict a nearly scale- 
invariant and gaussian distribution of fluctuations but 
higher order effects can break scale- invariance and gener- 
ate non-gaussianity |^|| . Defect perturbations are intrin- 
sically non-gaussian but are typically also scale-invariant 
in the generalized sense of "scaling" J?j]. We are primar- 
ily concerned here with the evolution of fluctuations from 
their initial state through the linear regime and do not 
consider these issues further. Note that changes in the 
spectrum of perturbations are simple to include in linear 
theory as evolutionary effects can be factored out into 
so-called "transfer functions" . 

B. Perturbation Type 

The perturbation type for the metric and matter fluc- 
tuations is the next most important distinction. A gen- 
eral linear fluctuation can be decomposed into scalar, vec- 
tor and tensor components. These manifest themselves as 
density, vorticity, and gravitational wave perturbations 
respectively and do not interact in linear theory. The 
scalar modes are the only ones that grow through grav- 
itational instability. Vector modes, on the other hand, 
always decay with the expansion. They can only be ac- 
tively generated by shearing (or anisotropic) stress in 
the manner. Tensor modes are intermediate. Left to 
themselves, they propagate as gravity waves, but they 
generate and can be generated by transverse-traceless 
(quadrupolar) stresses in the matter. 

The simplest inflationary models possess only scalar 
("S") fluctuations, and tensor fluctuations are generally 
cosmologically negligible in models with energy scales 



substantially below the Planck scale |8|]. However, mod- 
els whose initial conditions contain both scalar and tensor 
("ST") fluctuations are possible. Models with only S or 
ST generally have stresses that may be defined as func- 
tions of the metric, density and velocity perturbations 
and hence may be viewed as "passive" responses through 
equations of state. 

"Active" models have stresses that are a consequence 
of complex internal dynamics in the dark sector that 
cannot be simply specified as responses to gravitational 
perturbations. Although this definition is not precise 
for scalar and tensor modes, the very presence of vec- 
tor modes indicates an active source because these must 
be continuously generated to have an observable effect. 
Nevertheless, such models generally have all three types 
of perturbations ( "SVT" ) . 

C. Clustering Properties of Dark Matter 

Finally, the nature of the dark components affects the 
evolution of perturbations. We define as "dark" any com- 
ponent that interacts with the CMB photons only grav- 
itationally. Thus, even massless neutrinos are classified 
as dark matter in this scheme. 

Stresses in the dark components change the evolution 
of the mean density with time and the response of the 
matter to gravitational compression. We will loosely type 
models whose expansion rate is driven by a compress- 
ible type of matter (on scales relevant to cosmological 
structures) as "clustered" models and those which pos- 
sess matter that is incompressible as "smooth" models. 
We shall see that this distinction is in fact rather inex- 
act as it is not time invariant: essentially all models pass 
through phases when they would be considered smooth 
or clustered on the relevant scale. 

D. Phenomenology 

The key to understanding the phenomenology of a 
given model for structure formation is the evolution of 
metric fluctuations, in particular the Newtonian gravita- 
tional potential. Its qualitative behavior is determined by 
the initial conditions, perturbation type, and dark mat- 
ter content of the model. We illustrate this taxonomy 
scheme in Fig. |l|. 

The behavior of the gravitational potential is directly 
related to the evolution of density perturbations through 
the Poisson equation. Once its evolution is determined 
as a function of scale, not only is the present large-scale 
structure of the universe determined but also the whole 
time history of structure formation. The latter is im- 
portant for predicting the properties and abundances of 
high-redshift objects. 
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The gravitational potential also generates CMB 
anisotropies through gravitational redshifts |J and is the 
ultimate source of all anisotropies from scalar perturba- 
tions. CMB phenomenology can be essentially read off of 
the time evolution of the gravitational potential pi , al- 
though this involves understanding the backreaction from 
density perturbations in the CMB itself . 

Similar but simpler considerations apply for vector 
and tensor metric perturbations. They also generate 
anisotropies via gravitational redshifts but do not have 
unstable modes and hence do not affect large-scale struc- 
ture formation in linear theory. 

The difference between adiabatic and isocurvature 
models play a direct role in metric evolution because ini- 
tial curvature (or gravitational potential) perturbations 
are present in one and absent in the other. The per- 
turbation type changes the ratio of CMB anisotropies to 
large-scale structure. Finally, the dark matter proper- 
ties affect the evolution of the gravitational potentials. 
Smooth components by definition do not contribute to 
the gravitational potential but do contribute to the ex- 
pansion rate. They slow down the growth of structure 
and cause the gravitational potential to decay. Hence 
they decrease the amount of structure and increase the 
large-angle anisotropies of the CMB. 

In summary, the observable properties of structure for- 
mation models are encapsulated in the time evolution of 
the metric fluctuations. This in turn is governed by the 
stress properties of the matter both through its initial 
conditions and intrinsic properties. 

E. Current Model Zoo 

The archetypal model for structure formation is the 
standard cold dark matter model (sCDM), which is an 
adiabatic, passive, and clustered model. Here, an initial 
scale- invariant spectrum of adiabatic scalar ("S") per- 
turbations collapses via the gravitational instability of 
pressureless cold dark matter. Although this model is no 
longer viable from an observational standpoint, it pre- 
dicts phenomena sufficiently similar to the observations 
to act as a good starting point for model building. One 
of its failings is that it predicts too much small-scale 
power for the level of CMB anisotropies demanded by 
the COBE detection. 

A simple variation of the sCDM model that attempts 
to address this problem involves tilting the initial spec- 
trum of scalar perturbations (tCDM) to reduce small- 
scale power relative to large. Under certain inflationary 
scenarios, this brings about the addition of tensor per- 
turbations that further reduce small-scale density pertur- 
bations relative to the COBE detection. Such a model 
would be a "ST" variant of sCDM. 

The second class of variations involves changing the 
matter content so as to suppress the clustering of mat- 



ter, yielding a "smooth" variant. The prototypical ex- 
ample is the ACDM model, where an additional com- 
ponent of matter that does not cluster replaces most of 
the CDM. Another examples is the OCDM model where 
spatial curvature plays the role of the smooth compo- 
nent. Those two represent examples where the additional 
component is smooth on all scales and for all time by 
definition. Variants where the matter is only smooth 
on small scales include the HCDM model (e.g. |l2| also 
called C+HDM and MDM) with a component of hot 
dark matter, the (/>CDM model |l3| with a scalar field 
component 4> that tracks the background behavior of the 
matter, and QCDM JTJ,|l^] with a general scalar field 
("quintessence"). In a string-dominated universe (str- 
CDM) (l^] , the string network plays the role of a smooth 
component with the same equation of state as spatial cur- 
vature. Of course, one can have multiple smooth species 
as well, e.g. OACDM. The CDM class of models @ 
phenomenologically parameterizes all such models. 

Replacing the CDM with CDM of a different equation 
of state but no stress perturbations is a phenomenologi- 
cal possibility (GDM) suggested by |l7j] . This is an adia- 
batic, passive, and clustered variant of CDM. We will use 
the designation "CDMv" to represent all such variants of 
the CDM model. 

Isocurvature models have been proposed as alterna- 
tives to the CDMv class of models. The simplest ex- 
amples are those in which the initial stress fluctuations 
are established by balancing the density perturbations 
of two different types of matter. Examples include the 
axion isocurvature (AXI) [jl8| model and the primordial 
isocurvature baryon (PIB) model |l9| , where radiation 
density fluctuations are balanced by axions and baryons, 
respectively. The simplest versions involve only scalar 
fluctuations and hence are passive ("S") models. Models 
with and without smooth cosmological constant or spa- 
tial curvature components have been proposed. Versions 
with gaussian power-law initial conditions are observa- 
tionally challenged JfJ,|2l| but more complicated varia- 
tions exist po| . Based on our work, one of us has con- 
structed an isocurvature decaying dark matter (iDDM) 
model that defies conventional wisdom on isocurvature 
models and solves these observational problems [^|| 

Finally, topological defect models such as strings and 
textures fall into the isocurvature class but have fluctu- 
ations that are active ("SVT"). The simplest versions 
obey scaling and have only clustering matter but fail to 
generate enough large-scale structure for the observed 
CMB anisotropies |2^,^| . Models with a smooth A com- 
ponent have been proposed to alleviate these problems 
II- 

Hybrid models can also be constructed. If defects form 
after the inflationary epoch, one has a model with adia- 
batic initial conditions and active perturbations. A string 
model with inflation and cold dark matter (SIC) is a con- 
crete example p6|. One can also add in smooth compo- 
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FIG. 1. Taxonomy of structure formation. Models can 
be classified by their initial conditions (adiabatic or isocur- 
vature), perturbation type (passive or active), and cluster- 
ing properties of the dark matter (clustered or smooth on 
large-scale structure scales). Passive fluctuations involve 
stress responses to other perturbations and can support scalar 
("S") or scalar and tensor ("ST") components. Active stresses 
generate fluctuations and generally possess vector compo- 
nents as well ("SVT"). We will examine the extent to which 
this traditional categorization is useful in predicting model 
phenomenology. 

nents, e.g. spatial curvature (SICO). 

Clearly, the existing models do not even qualitatively 
exhaust the possibilities open to structure formation 
models. In the rest of the paper, we conduct an ex- 
amination of these possibilities beginning with general 
principles and explicitly stating the assumptions that are 
made in obtaining the models described as well as their 
generalizations. We summarize this analysis in a series 
of flowcharts (Figs. §, |, g, §, §, and ||.) 



III. COVARIANT PERTURBATION THEORY 
A. General Definitions 

We assume that the background is described by an 
FRW metric g^ v = a 2 ^^ with scale factor a(t) normal- 
ized to unity today and constant comoving curvature in 
the spatial metric jij . Here greek indices run from to 3 
while latin indices run over the spatial part of the metric: 
i,j = 1, 2, 3. The component corresponding to conformal 
time 



X° = 7 1 = 



dt 



(1) 



is 7oo = — 1 and 7oi = 7;o = 0. Unless otherwise speci- 
fied, overdots represent derivatives with respect to con- 
formal time and primes derivatives with respect to In a. 
c = 1 throughout. The background curvature is given 
by K = —Hq(1 — Sltot), where the Hubble constant is 
H = lOO/i km s" 1 Mpc" 1 . 

The ten degrees of freedom for the perturbations in the 
symmetric metric tensor g^ u can be parameterized as 

.9°° = -a 2 (l-2A), 



g — -a ±S , 



a 2 (Y 3 -2H Ll lJ -2H l S) 



T I ' 



(2) 



We refer to the lapse A as the potential, the three com- 
ponents of Bi as the metric shift, Hl as the curvature 
perturbation, and the five components of H l 2 as the met- 



ric shear following the conventions of 

Likewise, the symmetric stress-energy tensor can be 
parameterized by ten components 



J n 



-p-Sp, 
(p + p)(v l -B i ) 
+ P)v\ 
S P )S i j 



(P 



pU 1 



(3) 



i.e. the energy density and its perturbation (ft+ Sp), the 
isotropic stress (pressure) and its perturbation (p + 8p), 
the three components of the momentum density (p+p)vi, 
and the five components of the anisotropic stress tensor 
Hij. Note that the metric shift Bi enters in but not 
Tq 1 . Correspondingly, we shall see that Bi enters into the 
momentum but not the energy conservation equation. 

By writing the metric and stress energy tensor in this 
form, we have maintained general covariance. As a result, 
the equations of motion that result below take the same 
form for any coordinate system where linear perturbation 
theory holds. We reserve the term gauge invariant refers 
for objects that have the same value in each frame. 



B. Perturbation Representation 

While the perturbations are linear, they may be sepa- 
rated by their transformation properties under rotation 
without loss of generality (see Fig. ^). For covariant tech- 
niques that do not assume linear perturbations from the 
outset, see p7| , p8t and references therein. The five com- 
ponent metric shear and matter anisotropic stess 
n* 3 separate into one scalar, two vector, and two tensor 
components. The scalar stress/shear generates potential 
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FIG. 2. Scalar, vector, tensor decomposition. At the top 
of the tree of possibilities for structure formation models is 
the assumption that general relativity holds in the cosmolog- 
ical context and universe is homogeneous and isotropic in the 
mean with linear perturbations initially. Without further as- 
sumptions, the linear fluctuations may be expanded in scalar, 
vector, and tensor modes that do not interact while the fluc- 
tuations remain linear. 



flows (V x v = 0), whereas the vector stress/shear gen- 
erates vorticity (V ■ v — 0). The tensor stress generates 
tensor shear in the metric, which represents gravity waves 
in the transverse-traceless gauge. Therefore, the poten- 
tial A, curvature Hl, density Sp, and pressure Sp per- 
turbations are associated with scalar fluctuations alone, 
the metric shift Bi and velocity Vi with scalar and vector 
fluctuations, and the metric shear and anisotropic 
stress II^ with all three. Scalar, vector and tensor per- 
turbations may be treated independently in linear per- 
turbation theory. 

Fluctuations can be decomposed into the normal 
modes of the Laplacian operator Jy 



V 2 Q(°) = 




s, 




V 2 Qf 13 = 




v, 


(4) 


V 2 Qj 2 > = 




T, 



where vector and tensor modes satisfy a divergenceless 
and transverse-traceless condition respectively 



(±i) 



0, 7' 



)(±2) 



(±2) 



0. 



(5) 



In flat space, these correspond to plane waves times a 
local angular basis for the vectors and tensors 63] . 

Vector and tensor objects can of course be built out of 
scalar and vector normal modes through covariant differ- 
entiation and the metric tensor [HI 



(o) 



Q 



(0) 



(0) 



(6) 



)(±i) 



The perturbations in the kth eigenmode can now be writ- 
ten as 



A = iQ(°), H L = H L Q {0) , 
Sp = SpQW, Sp = 5pQ<®, 

which possess only scalar components, 



(7) 



Bi 



B^Q[ m) 

71= -1 



-(m) n {m) 
in — 1 

which possess scalar and vector components, and 

2 



(8) 



H Tij = £ 4 m) Q\? > 

m = -2 
2 

it, = ]T n< 



T (m) n (m) 



(9) 



m=-2 



which possess all three types. Here, scalar perturbations 
are denoted with superscript (0), which is elsewhere omit- 
ted. We will hereafter also omit the overhat in the normal 
mode amplitudes since real-space objects will no longer 
appear. Thus = v. 



C. Gauge Covariant Equations 

The equations of motion for the matter follow from 
the Einstein equations G ^ — QirGT^. Furthermore, the 
Bianchi identities guarantee TL, ;l/ — which represents 
covariant energy and momentum conservation. 

For the background, energy conservation implies 



P 



-3-(l+w)p, 
a 



(10) 



where w — pj p determines the background equation of 
state. Due to the isotropy of the background, momentum 
conservation yields no additional constraint. 

The Einstein equation determine the evolution of the 
scale factor through 



8ttG 



8ttG 



-a p c 



a 2 {p + Ps) ■ 



(11) 



Here we have divided contributions to the expansion rate 
into the ordinary density and an effective density com- 
ponent that does not participate in gravitational collapse 
and is hence labeled "S" for smooth. The curvature pro- 
vides the only component that is smooth by fiat 
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PS 



8irGa 2 



(12) 



with ws = —1/3. Even the cosmological constant is kept 
smooth simply by dynamics. However, this notation is 
convenient for considering components that are approxi- 
mately smooth. By keeping a general ps and p here, we 
avoid lengthy rederivation of the equations of motion for 
such cases. Note that p + ps = Per, the so-called critical 
density, and we can define a critical equation of state w cr 
by 



and the Einstein equations give 
{l-2K/k 2 ){kB {± ^ -H { T ±1] ) 



16TrGa z (p + p)(v (±1) -B (±L, )/k, (21) 
(kB^-H^) 



(±ih 



d „a 
drj a 



= -87rGa 2 pII (±1) 



(22) 



-3-(l + w C[ )p c 
a 



(13) 



Again the Bianchi identity reduces the number of inde- 
pendent equations to two. 

For the tensor modes, the Einstein equations reduce to 
a single relation 



Scalar matter perturbations obey the continuity and 
Euler equations 



drj 2 a drj ^ ^ ^ ^ 



H 



(±2) 



8irGa 2 pR 



(±2) 



d a 

;r + 3 - 

di] a 



3% 
a 



d 
drj 



a 



(P + P) 



6p 

(v - B) 



-(p + p)(kv + 3H L ), (14) 



(23) 



= fr-jj(i-3^)pn 

+(P + P)A, 



(15) 



and place 2 constraints on the 4 scalar matter- variables. 

The metric and matter are related by the Einstein 
equations 



(k 2 -3K)[H L + ~H T 



AitGa 2 



Sp- 



k 2 (A 



\h t ) 



= 87rGa 2 pn , 



±±(kB-H T )] 

3-(p + p)(v- B)/k 
a 

A + 2?) (kB H T ) 
drj a J 



1 



K 



A-H L - -H t - -^{kB - H T ) 



„ a n I a 
2- - 2 - 
a V a 



4irGa 2 (p + p)(v - B)/k, 

2 



a d 


k 2 ' 


A - 


d a 


a di] 


3 


drj a 



AnGa 2 (5p+ -5p) . 



(16) 



(17) 



(18) 



(19) 



Only two of these equations are functionally independent. 
The combination of these equations that corresponds to 
the conservation equation G^ v . v = is automatically sat- 
isfied by any choice of the 4 metric variables due to the 
Bianchi identities. The remaining degrees of freedom are 
related to gauge freedom as we shall see. 

Momentum conservation for vector perturbations gives 
the Euler equation 



d a 
— +4- 
drj a 



[(p + p)(u (±1) -B (±1) )/fc] 
-\{l-2K/k 2 )p\^ ±1 \ 



Neither the conservation equations nor the Bianchi iden- 
tity say anything about tensor perturbations. 

Although these relations are exact, they do not pro- 
vide a closed system. There are in general 10 equa- 
tions for 20 variables in the background and perturba- 
tions separately. For the background, homogeneity and 
isotropy brings this to 2 equations for the 3 variables 
(a), (p,p), where the grouping distinguishes metric and 
matter categories. For the perturbations, the general re- 
lations are broken up into 4 equations for the 8 variables 
(A, B, Hl, Ht), (Sp,Sp,(p + p)v,pH) for the scalar per- 
turbations, 2 equations for 4 variables (B^ 1 ^, H T ), 
(v^ 1 ), n^ 1 )) for each set of vector perturbations and 
1 equation for 2 variables (H T ±2 ^), (n( ±2 ') for each set 
of tensor perturbations. We can express the remaining 
1 + 10 degrees of freedom as the ability to choose the 
equation of state for the background w = p/p, the 6 
stress fluctuations (Sp^Tl^pTl^^pTl^ 2 ^) for the per- 
turbations, and the gauge (the 4 quantities (5rj, 6xi) for 
an arbitrary coordinate shift). 



D. Multicomponent Generalization 

The conservation equations (|l4|), ( |Ie| ) and ( ^p|) are 
valid for each species whose stress-energy tensor is in- 
dependently covariantly conserved. For example, they 
apply to the photon-baryon system and the dark sec- 
tor which only interact through gravity. The Einstein 
equations (|l^)-(|l9|), ( |2l| ) and ( |23| ) of course still hold 
with the appropriate summation over components, e.g. 
p = J2jP-J- Note that we do not include the smooth 
component ps in the multicomponent sum. 



IV. COORDINATE CHOICE 



(20) 
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A. Gauge Transformations 

The additional four component freedom in the Einstein 
equations is fixed by a choice of coordinates that relate 
the perturbations to the underlying smooth background. 
The most general coordinate transformation associated 
with the A;th normal mode is [|l| 



TQ 



(0) 



(24) 



T corresponds to a choice in time slicing and (L,L^\ 
L^ 1 )) a choice of spatial coordinates. Under the condi- 
tion that metric distances be invariant, they transform 
the metric as 



(1 



a 

w) — 

a 



_1 5pj_ 
PJ- 



(29) 



The common gauge choices of the next section all agree 
on the density perturbation in the clustered component 
well inside the horizon. 



B. Gauge Choice 

Gauge freedom can be used to simplify the equations 
of motion. Most commonly, it is employed to convert 
certain Einstein equations to algebraic relations and/or 
eliminate relativistic effects from the conservation equa- 
tions. 



A = A — T T , 

a 

B = B + L + kT, 

~ k a 
H L =H L --L--T, 
6 a 



Hn 



H 7 



kL. 



for the scalar perturbations and 

a(±i) = £<±i) + x,(±i) 



H, 



(±i) 



= H 



(±i) 



(25) 



(26) 



for the vector perturbations. 

Similarly, they transform the components of the stress- 
energy tensor as M 



Spj = Spj - pjT, 
dp j = Spj - pjT, 
vj = v.j + L, 



for the scalar perturbations and 



(±i) -(±1) 



r(±i) 



(27) 



(28) 



for the vector perturbations. All other quantities in the 
metric and matter are gauge invariant. In particular, 
the tensor modes do not exhibit gauge freedom since the 
transverse-traceless condition on Q^ ±2 ^ is sufficient to re- 
move the gauge ambiguity. The gauge is thus fixed by 
conditions on the metric which fully specify the transfor- 
mation (T, L, L^ 1 )) from an arbitrary frame. 

It is important to bear in mind that both the metric 
fluctuations (A, B, Ht, Ht) and the matter fluctuations 
(6p, Sp, [p + p]v) take on different numerical values in 
different frames even in this covariant notation. For ex- 
ample, if p evolves in time, a density perturbation 5p 
arises simply from the warping of the time hypersurface 
on which the perturbation is defined. Thus, a density 
perturbation differs negligibly only between frames sep- 
arated by (see Eq. JfO[) 



1. Vector Gauges 



Let us first dispose of the vector degrees of freedom. 



There are two natural choices flltl : H, 



(±i) _ 



= which fixes 



the gauge completely and B^ 1 " 1 = which leaves an 
arbitrary constant offset in H^ 1 ^. The latter does not 
produce a dynamical effect and can always be eliminated 
by specifying an initial condition for Hrf 



(±i) 



2. Comoving (Scalar) Gauge 

It is useful to consider a scalar gauge where the met- 
ric and matter fluctuations are simply related In- 
spection of the Euler and Einstein equations shows us 
that the coordinate choice B — v simplifies the equations 
of motion greatly. This fixes the time slicing through 
T — (v — B) jk. The additional condition Ht = speci- 
fies that L — HT/k and fixes the gauge completely. We 
call this the comoving gauge since here the momentum 
density vanishes. The remaining metric variables are la- 
beled A = £ and Hl = C This choice reduces the Euler 
equation to the algebraic relation for the potential 

{p + p)i = -5p+ 2 -(i-m/k 2 ) P n, (30) 

which is also simply related to the curvature through the 
Einstein equation (18) 



C= -Z + 4irGa 2 (p s +p s )v/k. 
a 



(31) 



Here we have again rewritten the curvature component 
as a smooth density contribution as in equation ( |TT| ) for 
easy generalization to approximately smooth cases. 

The simple relation between the metric and stress per- 
turbation of equations ( |30| ) and ( |3l| ) is what makes this 
gauge useful. A smooth component complicates these re- 
lations because of the difference between a frame that 



7 



is comoving with v versus the total-momentum-weighted 
velocity vp/(p + pg). With a constant comoving curva- 
ture, the continuity equation ( |l4| ) reduces to an ordinary 
conservation equation since metric changes to the fiducial 
volume are absent. 



As we shall see, employing both comoving and Newto- 
nian metric variables in the covariant language allows us 
to exploit the simple relations to comoving stresses in the 
former and comoving density perturbations in the latter. 



3. Newtonian (Scalar) Gauge 

Finally, the Newtonian gauge is defined by B = Ht = 
and labels A = \P and Hl = The gauge is completely 
specified through T = -B/k + H T /k 2 and L = -H T /k. 
The Einstein equations are reduced to algebraic relations 
that generalize the Poisson equation of Newtonian grav- 
ity 



(fc 2 - 3K)<S> = AirGa 2 



Sp + 3-(p + p)v/k 



fc 2 ($ + *) = -8nGa 2 pIl. 



(32) 



These algebraic relations and the fact that CMB 
anisotropics are simply related to $ and 4" make this 
gauge useful. 



4- Gauge- Covariant Variables 

It is often useful to speak of the variables of say the 
comoving gauge while in a Newtonian representation. 
Bai'dccn Q| introduced a so-called "gauge-invariant" lan- 
guage that achieves this. We denote such techniques as 
gauge covariant since they amount to introducing covari- 
ant expressions for objects that take on the desired mean- 
ing only in a specific frame. The only objects that cannot 
be made gauge covariant are those that are ill-defined due 
to coordinate ambiguities. 

To avoid confusion, we only use gauge-covariant vari- 
ables to describe metric fluctuations (C, £, <&, \&). Matter 
perturbations will always be represented in the comoving 
gauge unless otherwise specified. Note that v is the same 
in comoving and Newtonian gauges. 

The comoving curvature and density can be usefully 
expressed in Newtonian variables 



C = $ + 2 (* • 



*0*7 
P 



4irGa 2 5p= (k 2 - 3K)$ . 



(33) 



(34) 



obtained through equations ([25|) and (|18|). Likewise 
equation ( |3l| ) can be rewritten as 



Ps 



C - £ =(*-*') 
p 



(35) 



Recall that primes represent derivatives with respect to 
In a. 



V. STRESS PHENOMENOLOGY 
A. Stress Representation 

We have seen that the stresses of the matter compo- 
nents completely determine the evolution of perturba- 
tions (see Figs. ||, ||, and ||). The background stress is 
completely determined by the equation of state w. The 
scalar stress fluctuations are determined by functional 
relations between the pressure or isotropic stress pertur- 
bations Sp, anisotropic stress perturbation pIT and the 
density perturbation Sp. These relations may also in- 
volve hidden internal degrees of freedom. A model may 
also possess background stress without stress perturba- 
tions and vice versa. We call the former a smooth stress 
and the latter a seed stress. Vector and tensor pertur- 
bations likewise depend on the components pH^ 1 ^ and 
pUS ±2 ' of the anisotropic stress tensor. 

Scalar stress perturbations control the basic elements 
of the structure formation history, and so we pay partic- 
ular attention to categorizing their properties. As dis- 
cussed in §1VB, the isotropic scalar stress poses a spe- 



cial problem in that its value depends on the choice of 
coordinate or gauge. The comoving gauge (where the 
momentum density vanishes) provides a useful choice of 
gauge because of the simple relation between the metric 
and stress fluctuations. We will use these coordinates to 
define the total scalar stress 



S=-i = ^-\{l-3K/k 2 )- 
p + p 3 p 



-n. 



- p 



(36) 



It is useful to isolate gauge invariant aspects of the 
stress perturbation. The anisotropic scalar stress 



lirGa 2 pIl/k 2 



(37) 



is gauge invariant by definition; this form of the 
anisotropic stress also enters into the Einstein equations 
separately from S. 

An adiabatic stress perturbation obeys 



Sp A = (p + p)S A = —;5p. 

P 



(38) 



Although SpA is not gauge invariant, the adiabatic sound 
speed is 



c 2 = 



(39) 



Spa _ v[_ 
Sp p' ' 

This gauge invariance implies that there are no coordi- 
nate ambiguities when discussing the pressure support of 
adiabatic fluctuations. 
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The remaining gauge-invariant pressure perturbation 



Sr = —J— (tip- —bp 

p + p V p 



(40) 



Unlike adiabatic stresses, these may be present even when 
the comoving density perturbation dp is negligible (see 
Eq. Q]). Entropic stresses are the primary means of 
structure formation in most isocurvature models. 

The comoving-gauge analogue of Sa and Sr are also 
useful. As long as Sp/6p is less than unity, superhorizon 
stresses are negligible compared with curvature fluctua- 
tions. We call the part of 8p/Sp that is separable in time 
and space the comoving sound speed c 2 -, — f(k)g(r]) and 
the accompanying stress sonic 



Ss = c 2 c 



P + P 



(41) 



Separability is not a gauge-invariant property, but this 
is not in itself a problem because the comoving frame is 
dynamically special. 

We call the remaining isotropic stress the entropic 
stress 



S E = 



1 



P + P 

such that the total stress is 

2 k 2 



(6p - c 2 c 5p) , 



S 



S 



3K 



E 



3 8nGa 2 (p+p) 



Sn 



(42) 



(43) 



Note that if the comoving sound speed equals the adia- 
batic sound speed c 2 c — c 2 = p' j p' then Ss = Sa and 
Se = Sr- 



Scalar Perturbations 

Stress + Smooth — ► ACurvature 



Stress Free 

S, s n «C 



Scalar Stress 

S, 



FIG. 3. Scalar perturbations. It is useful to subdivide 
the stress-free class of scalar perturbations from the general 
possibilities. A "stress-free" perturbation has dimensionless 
stresses (S,Sn) that are much smaller than the comoving cur- 
vature perturbation £. Note that the stress-free perturbation 
condition does not preclude background or "smooth" stress. 



contribution that is truly external to the system of met- 
ric fluctuations. It is again possible to have large-scale 
entropic stress perturbation in the absence of initial cur- 
vature or density perturbation. 

The formal solution to equations (|44|) and (45) are 



(p s +Ps)v s /k = a -4 / ' dria A [5p a - |(1 - ZK / k 2 )p s Il s ] , 

Sp s = -a" 3 J d<na 3 [k(p s +p s )v s + 3^5p s ] . (46) 

The task of understanding a seed model like defects 
reduces to understanding its stresses, but this is a 
formidable task in realistic seed models such as cosmo- 
logical defects (e.g. p3fl). 



B. Seed Stress 

Seed stresses provide a special case with unique prop- 
erties. The effect of seed perturbations are in fact simpler 
to understand than those of the fluid type because the 
problem decouples completely. If the seeds do not inter- 
act directly with other types of matter, the conservation 
equations (14) and ( |l5| ) imply that the metric perturba- 
tions only affect the seed perturbations at second order 
since bare p s and p s terms may be dropped. They become 

M 



d „d 

;r + 3 - 

cirj a 



dp s = -k(p s + p s )v s - 3-Sp s , (44) 
a 



d , a 
+ 4 _ 

drj a 



{Ps+Ps)j = Sp s - H(l - 3K/k 2 ) Ps Tl s . (45) 



The basic principles of how stress fluctuations affect the 
gravitational potential still hold but here there is a stress 



C. From Stresses to Curvature 

These stresses are the fundamental sources and sinks 
of the metric perturbations. We therefore seek to express 
the comoving and Newtonian curvatures in terms of the 
stresses. Combining equations (p3|), (|35|), and (p6[) with 



yields 



and 



C' = -S + $ - S n ] 



Vp 



\t C + Sn. 
2 p 



(47) 



(48) 



(49) 



In a non-critical (non-flat, p ^ p cr ) universe, we combine 
these to eliminate Q from the evolution equation for the 
Newtonian curvature, giving 
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Vector Perturbations 




Stress Free 

H( ±i )»(p/p cr )n( ± D 



Anisotropic Stress 

H&D^Cp/p^nfii) 

V / 



1 



\ 



Pure Decay 

Decay from 
arbitrary initial 
conditions 



Stress Integral 

Defects 



FIG. 4. Vector perturbations. Vector perturbations simply 
decay from their initial value in the stress- free limit. The 
i ntegral solution in the presence of vector stress is given in 
§VTIIF and applies to defect models. 




Tensor Perturbations 




Stress Free 

H T (±2)»( P /p cr )ri(±2) 



Anisotropic Stress 

H T (±2)£(p/p cr )IK±2) 

V J 



\ 



\ 



Free Gravity 
Waves 

tCDM (matter) 



Stress Integral 

tCDM (radiation) 
Defects 



FIG. 5. Tensor perturbations. Tensor perturbations prop- 
agate as free gravity waves in the stress-free limit as is the 
case of a matter-dominated expansio n. The i ntegral solution 
in the presence of stresses is given in { VIII G and may be ap- 



plied to propagation during radiation domination as well as 
defect sources. 



EL 
p' 



2 p a 



2 Per 



S + Sji 



$' + 
1/4 

2 , 



P 



1^ 

2 Per 

p' 



_ El 
p' 

Sn ■ 



(50) 



Recall that the Newtonian curvature is simply related 
to the comoving density perturbation through equation 



In a critical-density universe (p = p cr ), equations (|46 
and (|49[) can be formally solved as integrals over the 
stress fluctuations to yield 



C(a) - C(0) - 



da 



S, 



(51) 



and 



* = C-^ f^[C-S~S n]+ C^. (52) 
a J y/p a 

The last term is the decaying mode of $ where C =const. 

There are three general conclusions that we can draw 
from equations (pj) and (M). The first is that in the 
absence of an initial comoving curvature perturbation 
(C(0) = 0), a stress fluctuation will generate one of or- 
der £ ~ —S. The same goes for the Newtonian cur- 
vature $ ~ —5. The reason for this behavior is that 
a stress gradient kSp generates a potential flow with 
(p+p)v ~ (krj)5p which generates a density perturbation 
of Sp ~ —(kri) 2 8p and hence a curvature perturbation of 
$ ~ —5p/p cr . Note this intuitive argument fails for other 
gauge choices. 

Second, starting with a curvature perturbation and as- 
suming sonic stresses Sp = c^(5p, it is clear that the same 
mechanism of generating flows will generate an oppos- 



ing curvature fluctuation A<f> 



7 5p/p 



that will destroy the initial curvature fluctuation when 
cckrj ~ 1. In physical terms this occurs because pres- 
sure support prevents perturbations from collapsing and 
hence causes the curvature perturbation to redshift away. 

Finally, the anisotropic stress contributes to the to- 
tal stress S and thus can both create and destroy co- 
moving curvature fluctuations. Furthermore, it enters 
separately into the Newtonian curvature through equa- 
tion (p2|) . This is because the Newtonian frame, unlike 
the comoving frame, is defined to be globally shear free 
(B = Ht = 0). The coordinate transformation that 
maps the comoving frame to the shear free frame depends 
on the anisotropic stress and hence aliases background 
evolution into contributions to the Newtonian curvature. 

Unfortunately, despite their general appearance, equa- 
tions (|5l]) and (|52|) are only formal solutions since the 
time evolution of the stress sources generally depend on 
the curvature fluctuation itself. We will use the spe- 
cial properties of smooth, anisotropic, entropic, and sonic 



stresses to address this problem in §VI, § VII, and § VIII 



D. From Curvature to Observables 



As discussed in SIID 



the Newtonian curvature is di- 
rectly related to observables in the CMB and large scale 
structure. We are now in a position to quantify these re- 
lations. The Newtonian curvature $ and potential ^ en- 
capsulate all observable properties of scalar fluctuations. 
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The contribution of a given £;-mode to the amplitude of 
£th multipole moment of the CMB anisotropy is given by 



21- 



- w P drie- T \ [* - $ + f(9 
1 Jo L 



<je[k(vo - ??)] + Tv b j' e [k(r] a - 77)]} , (53) 



where we have dropped the small correction due to the 
polarization of the CMB. r is the optical depth to Comp- 
ton scattering between the present (779) and the epoch in 
question (77); O is the photon temperature perturba- 
tion in Newtonian gauge, and Vb is the baryon velocity 
in Newtonian or comoving gauge. For an open universe, 
the spherical Bessel function jj> is replaced by the hyper- 
spherical bessel function. Note that the prime here and 
here only refers to a derivative with respect to the ar- 
gument of the Bessel function. With the random phase 
assumption for the fc-modes, the scalar contribution to 
power spectrum of the anisotropies is 



C, = 



dk 



(2£+l) 2 



(54) 



In equation (J53J) , the ^ — $ term leads to the so-called 
integrated Sachs- Wolfe (ISW) effect and contributes once 
the optical depth to scattering becomes small, i.e. after 
last scattering. The other terms, Go + 'J' and v bl are local- 
ized to the last scattering surface itself and represent the 
effective temperature of the distribution and the Doppler 
effect respectively. They are responsible for the so-called 
acoustic peaks in the CMB spectrum, the morphology of 
which can be directly read off of the metric driving terms 

The influence of time variability in and $ depends 
on how the variation rate 1/A?7 compares with the per- 
turbation crossing time for sound (before last scattering) 
and light (after last scattering). For variations on a much 
shorter time scale (kAri < 1 or kc s Ar/ <C 1), only the 
total change A(4" — <£) is observable since all photons 
suffer a uniform gravitational redshift. For variations on 
a much longer timescale the effects mainly cancel out as 
the photons either traverse many wavelengths of the fluc- 
tuation during the variation with redshifts and blueshifts 
cancelling or equivalently undergo many acoustic oscilla- 
tions. Changes whose duration are synchronized with the 
oscillation period are the most effective. These general 
considerations apply to metric variations of the vector 
and tensor type as well. Two commonly encountered ex- 
amples are the cancellation cut off at high £ for a uniform 
potential decay and the driving of acoustic oscillations 
from synchronized potential decay in the radiation dom- 
inated epoch. 

A constant potential also leads to observable effects 
through the effective temperature term Go + "P. Since 
Oo is the temperature perturbation in Newtonian gauge, 



it can be obtained by a gauge transformation from the 
comoving gauge by noting that in that gauge both the 
density perturbation and the potential £ are negligible 
because stress perturbations must be negligible for the 
potential to remain constant. Equation ( |25| ) and (|27 
then imply 



S P-r " /; 
6 = — - - -v/k, 

Ap^ a 

* = v/k + -v/k . 
a 



(55) 



For adiabatic fluctuations, <5/? 7 /p 7 oc Sp/p and is hence 
negligible outside the horizon by virtue of the Poisson 
equation ([54]). Since $ is a constant by assumption, these 
equations can be integrated to give 



6 



3 1 



(56) 



and hence 6 + * = *(1 + 3w)/(3w + 3). This rela- 
tion ultimately comes from the fact that ^> represents a 
time shift and a oc £ 2 /3(i+™) |rj. For w = 0, this re- 
duces to the well-known result that the effective temper- 
ature is ^/3 in the adiabatic sCDM model ||. Compared 
with this model, those that are dominated by the ISW 
term — $) w 2^ — Su like traditional isocurvature 
and smooth models potentially have up to 6 times the 
anisotropy for a given potential fluctuation. 

The behavior of the density perturbations that under- 
lies large-scale structure are even more directly related to 
the Newtonian curvature perturbation. Inside the hori- 
zon, all reasonable choices of gauge agree on the density 
perturbation. In particular, the comoving gauge den- 
sity perturbation is algebraically related to $ by equa- 
tion ( |34| ) and hence allows a simple translation of results 
for one to the other. Consequently, the relation between 
temperature and potential fluctuations discussed in the 
last paragraph translates directly into a relation between 
density perturbations and CMB anisotropies. 

With this relation, the so-called transfer function of 
the density perturbations below the current horizon can 
be simply read off of the time history in the potential 
without the usual gauge ambiguity in defining the initial 
density perturbation, 



T(k) 



<S>(yo,k) $(0,0) 
$(0,fc) $(770,0) 



(57) 



The power spectrum of density perturbations today is 
then proportional to T(fc) 2 Pj n i t i a i- Any process that 
makes the potential decay relative to the k = mode 
will produce a downturn in T(k) and a reduction of small- 
scale power. 
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OCDM 
HCDM 



(|>QCDM 
AHCDM 



strCDM OHCDM 
iSfCDM QHCDM 
QCDM 
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FIG. 6. Stress-free scalar perturbations. II all ol the matter 
which drives the expansion participates in gravitational in- 
stability, the matter is said to be "clustered" otherwise there 
exists a "smooth" component. In either case, the perturba- 
tions depend only on the background ("backgnd.") equation 
ol state w and can be tracked with simple techniques (mid- 
dle row), which lead to exact solutions (bottom row) that 
describe behavior in a wide range ol models. 



VI. STRESS-FREE PERTURBATIONS 

Adiabatic models for structure formation, those whose 
initial conditions contain true comoving curvature per- 
turbations, all go through a period in which the scalar 
stress perturbation may be neglected. This is a direct 
consequence of causality. Stress gradients affect struc- 
ture formation simply by the causal motion of matter. 
On scales larger than the horizon, the change in the co- 
moving curvature perturbations due to causal motion al- 
ways can always be neglected. 

We begin in §VIA with adiabatic models that have 



no smooth components. Although the comoving curva- 
ture remains at the value set by the initial condition, the 
Newtonian curvature depends on the equation of state. 



We consider then in § VI B the effect of a smooth compo- 



perturbation). For each case, we begin with a general 
description of the resultant phenomenology. We then 
illustrate the phenomenology with full solutions of the 
perturbation equations and discuss applications within 
the current generation of structure formation models of 
§ II E . We will follow this pattern throughout the paper. 



An overview of results is given in Fig. 



A. Clustered Case 

We begin with the case in which there is no smooth 
component (ps = 0) and stress fluctuations are negligi- 
ble compared with metric fluctuations (S <C C)- Here, 
equation (^8|) simply implies 



( = const. 



(58) 



Equation (p2|) then gives 



4> 



,[l-&[±) + &[*±S a + C&. (59) 
a J VP/ a J VP « 



We will focus only on the behavior of the first term. 
The second term may be neglected if Sri <C C- However, 
unlike S <C Sn *C C is n °t a consequence of causality 
and is mildly violated in the case of free radiation |Tc| ] 
and can be strongly violated in defect models [|lj. We 
return to consider its effects in ^VII A . The last term 
is a decaying mode that carries no comoving curvature 
(C = 0). For a constant equation of state w, it scales as 



$ OC o -3(l-Mi0/2- 



(60) 



It is apparent from the form of the first term that since 
p cannot grow with time (w > —1), the first integral in 
equation (|59j) goes to a constant between and 1 with 
the two extremes representing w — > oo and w = — 1 re- 
spectively. The integral is dominated by the most recent 
epoch; only the equation of state at the epoch in question 
matters. It is simple to show that during periods where 
w is approximately constant ]32f| , 



3 + 3tu 
5 + 3w 



(61) 



At each epoch where w decreases from w± to W2, the 
Newtonian curvature decreases by 



$(wi) - $(w 2 ) _ (wi - w 2 ) 

$(t«i) ~ (1 + wi)(5 + 3w 2 ) 



(62) 



nent on the Newtonian curvature (and hence the density 



These results are easy to understand in the non- 
relativistic limit. In the absence of stresses, the gradients 
in the gravitational potential set up flows as w ~ (kr])$. 
The divergence of this flow generates density perturba- 
tions S ~ (1 + w)(krj) 2 § for constant w. By the Poisson 
equation $ ~ (krj)~ 2 S, this is exactly the rate of growth 
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needed to keep the potential constant. Here we have used 
the fact that A-kGo 2 p = 3(a/a) 2 /2 ~ r\~ 2 if there are no 
smooth components. When w decreases from w\ to W2, 
the change affects the proportionality constants but not 
the scaling between 5 and <£>. This leads to a decrease in 
the Newtonian potential to a new constant. 

Full Solutions. — If the matter is predominantly com- 
posed of two components with different but constant 
equations of state w\ and 11)2, we can solve equation ( pHT 
exactly provided Sn<(- The result is 



3wi 



-F 



y 



l + y 



(63) 



= 1 - 



where n = (5 + 3wi)/(6u>i — 6W2) and 
y = P2 /pioca 3 ^- w ^. 



(64) 



Here, B x (p,q) is the incomplete beta function and 
F(a, 6; c; x) is Gauss's hypergeometric function (i.e. 2-Pi)- 
The combination y/ (l+y) often enters into such solutions 
as it is just P2I Per-, the fractional density perturbation 
supplied by component 2 as a function of time. 

Any case with 2n equal to an integer can be expressed 
in elementary form. In particular, ifn = ./Vor.ZV+l/2 
for an integer N > 0, then 



1; x 



(65) 



1 T{n + l)T(n- k- 1/2) 



k=0 
X - 



Y(n - fe)r(n + l/2) 



x - 1 



N 



T(n + l)T{n- N +1/2) 
T(n- N +l)T(n+ 1/2) 



fn 



where 




1 log^, 



N; 

N - 



(66) 



Applications. — These solutions apply to any adia- 
batic model with a matter-dominated epoch and may be 
used to explore the behavior of perturbations entering 
and exiting the matter-dominated epoch. For example, 
the matter-radiation case reduces to [p3| 



$ _ 3 2 8 16 
~( ~~ 5 + 1% ~ 15y 2 ~ 15y 3 



1671+^ 
15y 3 



(67) 



This solution only holds in the absence of anisotropic 
stress perturbations and does not strictly apply to 
the usual radiation components. The neutrinos carry 
anisotropic stress [lOj as do the photons after recombi- 
nation. We consider how such effects can be taken into 



account in SVIIA 




0.2 0.4 0.6 

y/(i+y) 

FIG. 7. Stress-free clustered case. In the case where all 
contributions to the expansion rate also cluster, the comoving 
curvature is constant but the Newtonian curvature changes 
with the equation of state w. For the case of the transition 
between wi — and W2 between —1 and —1/6, we show here 
the analytic solution of equation (^) is compared with a full 
numerical Boltzmann solution of a QCDM model (including 
radiation) for a mode that is outside the horizon at the given 
epoch. Notice that only for W2 = — 1 (cosmological constant) 
does the potential decay to zero as p%l 'p CI = y/(l + y) — > 1. 

For adiabatic models where there exists a component 
with W2 < |l4| , |i5| , |l7t , these solutions describe the exit 
from the matter-dominated epoch. A comparison of the 
analytic solution with full numerical solutions in those 
cosmologies is given in Fig. [?]. A special case is w± — 
and u>2 = —1, the matter to cosmological constant 
transition, where Ip5l,p6i 



1 - \(l + y) 1/2 y- 5/e B y/(1+y) (5/6, -1/3) 



(68) 



As is evident from equation (|61j), the cosmological con- 
stant is the only case where the Newtonian curvature de- 
cays to zero. Note in particular that W2 = — 1/3 does not 
correspond to the behavior of spatial curvature in spite 
of the fact that the effect on the expansion rate is the 
same. The presence of fluctuations in the W2 component 
prevent the gravitational potential from decaying to zero. 
Across the W\ = (matter-dominated) to W2 = —1/3 
transition, the Newtonian curvature goes from 3£/5 to 
C/2. Thus, the string-dominated (strCDM) model which 
has w s tr = — 1/3 does not behave like an open model on 
the large scales relevant for the CMB (c.f. Jig]) 



B. Smooth Components 

We next consider how smooth components affect the 
growth of structure. Although curvature is the only com- 
ponent that is smooth by definition and a cosmological 
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constant the only one that is smooth by dynamics, under 
certain circumstances other components can be approxi- 
mately smooth. 

We define a component (S) to be smooth if its den- 
sity fluctuations are small in comparison to those of 
the clustering components (C), i.e. Sps <C Spc re- 
gardless of whether ps < pc- To be maintained dy- 
namically, the respective energy fluxes must also sat- 
isfy (ps + Ps)vs <C (pc + Pc) v c- ln this section, we 
add the subscript (C) on the remaining matter to re- 
mind the reader that part of the total matter density has 
been designated as effectively smooth in the division of 
equation (|ll|). Our analysis also applies to cases where 
dps <C Spc where the time-ave rage is over the dynamical 
time of the C component (see VIII C ). 

One cannot demand that Sps be identically zero since 
the continuity equation ( |l4[ ) generates a density fluc- 
tuation as the metric curvature changes unless ws = 
Ps/ps = — 1 ( a cosmological constant). The fractional 
density fluctuation is generically at least of order the cur- 
vature fluctuation Since this term exceeds the energy 
flux term outside the horizon due to causality, energy 
conservation forbids smooth components (with ws =/= 1) 
on these scales. 

Components that are smooth within the horizon are 
possible. For example, the C component may be driven 
to collapse by potential gradients while the S compo- 
nent is supported against collapse by stress gradients 
(see in §VIIIC). Here Sps/ps ~ $ *C Spc/pc- Since 



these stress gradients are set up exactly so as to keep 
the component smooth, one can replace such stress ef- 
fects in equation ( |50| ) with an additional smooth density 
component p4| . The remaining perturbations can then 
be approximated as stress-free and generate curvature 
fluctuations as 



1 



P'c 



2p c 



2 



Pc 



0. 



(69) 



Now even for constant wc, = no longer solves 
the equation of motion. Mathematically, ps adds to 
the expansion drag (<!>') terms but not the gravita- 
tional (<E>) terms. Physically, potential flows still cre- 
ate density perturbations as pc ~ (pc + Pc)(ki]) 2 Q 
but the Poisson equation leads to a smaller potential 
$ = (kr])~ 2 Spc/(pc + Ps)- As this process continues, 
$ decays away. 

Equation ( |69| ) has simple solutions in the limit that 
p' cx 3> p' c , as is usually the case when the smooth com- 
ponent dominates the expansion. In this case, the general 
solution to the equation is 



$ = C x a 1 + C 2 a 



din a 



aPc 

1/2 
Per 



(70) 



Both terms represent decaying modes as long as w a — 
2w c < 1. 

Full Solutions. — The full solution to equation ( |69| ) can 
be obtained analytically for a clustering component with 
a constant equation of state and present-day fractional 
density contribution pair (wi,Clci), a smooth component 
with (wi,Qsi), and/or a second smooth component with 
(w 2 ,^S2)- 



*j = y 



_ „l/3Atu [ v_ 



1+y 



(71) 
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fisi + &C1 



(72) 



where j = 1,2 for the growing and decaying modes corre- 
sponding to — , + in the ctj equation. Note that negative 
contributions to the critical density, e.g. positive spatial 
curvature K > 0, are also covered by these solutions. 

The hypergeometric solution allows one to identify ele- 
mentary solutions more easily. In particular, cases where 
ay = — j7Vi/2|, (3tui + l)/6Au> = N u or 2a, - (3wi + 
l)/6Aw + 1/2 = N% can be expressed in terms of ele- 
mentary functions. Cases of the form F(a, b; b + Ni; x), 
F(a,b;a + Ni;x), or F(N x /2, N 2 /2; N 3 /2; x) can be ex- 
pressed in terms of elementary functions, incomplete beta 
functions, and/or complete elliptic integrals. Here, the 
Nj are integers and a and b are real numbers. For exam- 
ple, the case of w% — and w% = —1/6 can be expressed 
in closed form for any Oci/C^Ci + ^si) : 



oc 



y a > 2^ + i + yr+y 

a (l + VTTy) 2 ^ +1 



(73) 



For w\ = 0, any case in which 2a j + l/6u>2 + 1/2 is 
an integer can be simplified. One can also simplify the 
growing mode of cases with w\ = flsi = and w^ 1 equal 
to an integer. 

Finally for completeness, there is a well-known special 
case of equation ( |69| ) that is not completely covered by 
equation ([FT]) but that does have an integral solution. 
This involves a clustering component w — (CDM) and 
a smooth component composed of an arbitrary admixture 
of w = — 1/3 (curvature) and w = — 1 (A) pieces p7||: 



$ a a p, 



- 1 n 1 / 2 



din a 

9 3/2' 
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(74) 



for the growing and decaying modes respectively. Of 
course, if either the curvature or A component is neg- 
ligible, the solutions have an analytic form described by 
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0.4 0.6 

y/(i+y) p 

FIG. 8. Single smooth component. The same as Fig. M 
except that W2 is taken to be a smooth component. Here 
equation ( fri| ) is compared to a QCDM model but for a mode 
that is well inside the horizon at y = 1. Discrepancies at 
early times are due to radiation contributions in the QCDM 
models and are particularly pronounced for W2 — > since then 
the epoch when the Q and matter components are equal is 
driven into the radiation dominated era for this flat Q. m = 0.35 
model. 



equation (|71|). Unfortunately, this solution cannot be 
generalized to arbitrary combinations of smooth compo- 
nents since it relies on the fact that w = —1/3 does not 
accelerate the expansion and that w = — 1 gives a con- 
stant density contribution. 

Applications. — Smooth components are widely found 
in adiabatic models and their description in terms of 
equation ([n]) are given in Table |. 

In fact, all CDM variants (CDMv) go through a phase 
in the radiation-dominated epoch when the perturbations 
in the radiation are pressure supported leading to an es- 
sentially smooth component S7 rac j = £ls2 with u>2 = 1/3 
and also a smooth baryonic component Of, = f2gi with 
wi = that is held against collapse by Compton coupling 
to the photons (see §VIIIC and J3l| ), In this case, the 
clustered matter is the CDM (JIcdm = ^ci)> an d the 
solution ( fn| ) describes the evolution of the CDM density 
perturbations via the Poisson equation. The behavior 



deep in the radiation domain is given by equation (70), 
which says that density perturbations actually grow log- 
arithmically. The full solution ( fn) ) maps this logarith- 
mic mode into a power-law growth across the matter- 
radiation transition and is useful for determining the am- 
plitude of small-scale fluctuations as a function of the 
baryon content Q . 

Smooth components that dominate at late times are 
also described by equation ((7^) . Many such models have 
been proposed to reduce the amount of small-scale power 
in the standard CDM model. The prototypical exam- 



ples are the ACDM model where f2g2 = = 1 — 

(1V2 = — 1) and the open model OCDM fls2 = 

1 — Q m (w2 = —1/3). The former case is special in that 

it may alternately be considered a clustered component 

(see § VI Ap . The latter case must be considered as a 

smooth component, and the growing mode reduces to 

M 
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(75) 



where y oc a. 

The obvious generalization of such models involves a 
smooth component with an equation of state that dif- 
fers from the curvature or cosmological constant exam- 
ples. The prototypical case is the HCDM model where a 
massive neutrino (or hot dark matter) component with 
101 K remains smooth on small scales due to residual 
relativistic effects. Our general solution in fact allows an 
additional smooth component W2, which could be cur- 
vature (OHCDM), a cosmological constant (AHCDM), 
or even some new component with an equation of state 
w 2 = wq (QHCDM). In fact, in the case where there 
is no hot component, the latter solution describes the 
"quintessence" model (QCDM) that has recently received 
much attention Jl4],[l5]]. Here, a scalar field supplies a 
density component that is smooth inside a sound horizon 
that corresponds to the particle horizon. A comparison 
of the analytic solution to numerical results for QCDM 
and AHCDM are given in Fig. ^ and |^. Likewise, the so- 
lution applies to the GDM generalization |l7j] where the 
sound horizon is allowed to be arbitrary. 

Finally, equation ( |74| ) describes the case of a CDM 
model with both curvature and cosmological constant 
contributions. 
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TABLE I. Correspondence between the analytic solution 
of equation ( |7l| ) and models with smooth components. Q.k 
is the fractional effective density supplied by the curvature 
component. 
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FIG. 9. Two smooth components. Equation ( |7l| ) is com- 
pared with numerical solutions in a AHCDM which has 
smooth A and hot (H) dark matter Q v components. Nu- 
merical solutions are for modes well below the Jeans scale of 
the hot dark matter and discrepancies at early times are due 
to radiation contributions in the model. 




X = k?7 

FIG. 10. Gravity wave modes. The amplitude of free grav- 
ity waves remains constant outside the horizon and oscillates 
and decays inside in a manner dependent on the equation 
of state of the background. Shown here is the derivative of 
the mode since it acts as the source of radiation anisotropies. 
Note that as w increases the oscillatory phase begins sooner 
relative to horizon crossing. 



C. Vector Perturbations 

The behavior of vector modes is far simpler than that 
of scalar modes. Vector anisotropic stress can be ne- 
glected if Hj^ 1 '* 3> pll^ 1 ) I p in a gauge where B^ 1 " 1 = 0. 
In the absence of stress perturbations, vector modes sim- 
ply decay with the expansion (see equation [gfj ) is solved 
by 



(p + p)(w (il) -B (±1) )oca- 4 



(76) 



The metric source is related algebraically to this quantity 
by equation (|2l|). 



D. Tensor Perturbations 

( ±2) 

The tensor anisotropic stress is neglig ible if H^> » 
pIT ±2 ) / p CI . In this limit, tensor metric fluctuations re- 
main constant outside the horizon regardless of the ex- 
pansion rate and propagate as free gravity waves inside of 
it. They are described by equation (^3|) with the sources 
set to zero. 

Full Solutions. — If the expansion is dominated by a 
component with constant equation of state w CT > —1/3, 
the fundamental modes of gravity waves are 

gl = 2 " + ' r ^ +3/2) ^„M. 

a. g^tffl ,^, (77) 



with x — kr), m = (1 — 3w)/(l + 3w) and fc = \lk 2 — 2K. 
Note we have normalized the modes so that -Hi(O) = 1; 
this mode drops from unity into damped oscillations 
when the wavelength reaches some fraction of the hori- 
zon that decreases with m and hence increases with w. In 
Fig. [To], we plot the derivative of these modes, as that de- 
termines its effect on the radiation through gravitational 
redshifts. 

For w < —1/3, the universe accelerates and the hori- 
zon stops growing with the scale factor. This implies 
that gravity waves will also freeze out at some finite value 
related to their amplitude when the universe began ac- 
celerating. The solutions for the gravity wave behav- 
ior relative to the epoch of freeze-out follow the form of 



equation (77) with x = k dij. The Hi mode then has 
damped oscillations as x increases from negative values 
but freezes in to finite value as x — > in the infinite 
future. 

Applications. — The solutions above help us under- 
stand the phenomenology of tensor anisotropies in the 
CMB. The tensor analogue to equation (|5^ ) is 




(78) 



where A77 = rjo — r\. The power spectrum is again given 
by equation (p4|). The appearance of the e~ T damping 
reflects the fact that anisotropic stress cannot be sup- 
ported in the optically thick limit; we take r — > in 
the examples here to consistently neglect all anisotropic 
stress effects. This in fact is a good approximation for 
tensor anisotropies in the neutrino background radiation. 
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FIG. 11. Gravity wave anisotropies generated in free radi- 
ation. Numerical solutions are from a Boltzmann code with 
scattering sources removed. Notice that as w increases so does 
the anisotropies. This is directly related to the behavior of 
the gravity wave modes in Fig. hoi 



For a flat universe with constant w > —1/3, the results 
are shown in Fig. [ll] for the same scale-invariant initial 



spectrum of gravity waves, k \H, 



(±2) 12 



=const. Notice 



that the anisotropies decrease as w is decreased. Like the 
scalar ISW effect discussed in W D, the contribution of a 



decaying tensor mode to the anisotropy depends on how 
long the gravity wave takes to decay relative to the light 
travel time across the perturbation. In the limit that it 
decays before horizon crossing, the photons experience 
the full gravitational effect. If it decays well after hori- 
zon crossing, then the effect suffers cancellation as the 
photons traverse many wavelengths of the perturbation. 
Thus, the relative contribution to the anisotropy can be 
read off the behavior of the normal mode in Fig. [l(]. 
As w increases, changes in Ht occur at smaller times 
relative to horizon crossing. The anisotropy contribu- 
tion accordingly goes up. The effect is most dramatic 
for the quadrupole since all A:-modes above the horizon 
contribute to the quadrupole. This effect explains why 
the tensor spectrum in the usual matter-radiation uni- 
verse shows an upturn in the spectrum as one goes from 
modes that crossed the horizon in the matter dominated 
epoch to those that crossed in the radiation dominated 
epoch (see Fig. |l7| ). 

For w < —1/3 the gravity waves freeze out. Since 
CMB anisotropies are driven by changes in the gravity- 
wave amplitude, the addition of a w < —1/3 component 
should suppress anisotropies; this prediction is in agree- 
ment with the effect found in ACDM and QCDM models 



VII. PURE STRESSES 

We have shown in the previous section that in certain 
regimes stress perturbations can be ignored. However, 
to provide a complete history of structure formation, one 
must track the perturbations across all scales and time. 
We will first consider the pure stress cases in which the 
dominant stress contribution is anisotropic, entropic, or 
sonic, as defined in 



VA 



These prototypical cases have 
analytic solutions and are the starting point for the gen- 
eral cases discussed in § VIII . They also have direct appli- 



cation in many models. The anisotropic and sonic solu- 
tions are applicable to all adiabatic variants of the CDM 
model (CDMv). The entropic solutions show how pro- 
totypical isocurvature models such as the baryon (PIB) 
or axion (AXI) isocurvature models form structure. We 
outline these results in Fig. 12. 



A. Anisotropic Stress 

Anisotropic stresses play a special role because they 
enter directly into the Newtonian metric through Su, as 
opposed to other stresses which only contribute through 
the causal motion of the matter. This is despite the fact 
that the comoving curvature £ depends only on S. In an 
isocurvature model, only the comoving curvature need 
vanish initially, not the Newtonian curvature. 

Similarly, even though £ remains constant above the 
horizon in adiabatic models, the Newtonian curvature $ 
evolves under anisotropic stresses (see Eq. j52|). Recall 
that 



Su 



iirGa 2 pIl/k 2 



-8|SV^n/*», 

a J Per 



(79) 



such that its effect is enhanced by [p 1 / p cr )(kr])~ 2 rela- 
tive to S. Note, however, that once the universe en- 
ters a period when (p/p CI )H <C (ki]) 2 (, all traces of the 
anisotropic stress from any previous period vanish in the 
relation between $ and £. Its effect does not vanish 
from the CMB, however, since anisotropies record a time- 
integrated history of the gravitational potentials. 

Full Solutions. — To close this system of equations, we 
need a relation between Su and £. We will consider two 
limiting cases: when Su < C as i n the case of stresses from 
radiation backgrounds, and when Su S> C as is possible 
with models involving active sources such as defects. 

In the former case, anisotropic stress is generally cre- 
ated as a by-product of gravitational instability. Its 
anisotropic nature suggests that it can be created from 
shear in the velocity and metric, and its coordinate trans- 
formation properties demands that its source be gauge 
invariant. The linear combination of these sources that 
satisfies these requirements is (kv — HT)- The anisotropic 
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FIG. 12. Scalar stress perturbations. Scalar stress perturbations are divided into three pure classes: anisotropic, entropic, 
and sonic. If one of these is dominant, then general techniques (third row) can be applied to analyze the resultant behavior of 
fluctuations. Many models (fourth row) do go through phases where the stresses are dominated by one of the pure stresses. 
The designations "all models" and "all adiabatic models" assume a neutrino background and fluctuations that were present 
before last scattering. Several cases with analytic solutions illustrate the range of behaviors. 



stress can also have a dissipation timescale rju- Together 
these considerations imply an evolution of the form 



ri + ? /n 1 n = 4(kv - H T )a. 



(80) 



In this section, we are interested in large-scale effects and 
hence want the longest timescale for the dissipation; this 
is set by the expansion time. We will consider cases where 
the timescale is much smaller than the expansion time in 



VIII D 



If we take r] u — 3a/a, we recover the phenomenologi- 
cal parameterization of [p7|p| 



n + 3-n = A(kv - H T )a. 
a 

which has the formal solution 



n = 4a" 3 J dria 3 (kv - H T )a . 



(81) 



(82) 



*a = Cvi s /w 9 in the notation of 117 



From equation 



we find 



i J 1 + w p \a 



a 
a 



(83) 



for modes well under the curvature scale. Employing 
equation (61) for the zeroth-order potentials and assum- 
ing constant oc and p = p cr yields 



7 = ~ 2a{kll) (4 + 3«,)(5 + 3«,) ' 



(84) 



for constant w. With this source in equation (|52j), the 
Newtonian potential becomes 



$ 3 + 3w 



C 5 + 3w 
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FIG. 13. Anisotropic stress at large scales. The presence ol 
anisotropic stress affects the Newtonian ($) but not comoving 
(£) curvature. II Sn < C solutions for $ can be constructed 
by iteration. Here we show an example where the anisotropic 
stress is produced by massless neutrinos and compare numeri- 
cal results from a Boltzmann code with the iterative solution. 
To first order, Sn/C = 4fV„/15(4.4 + N„), where N v is the 
number of massless neutrino species. 



if the curvature contributes negligibly to the expansion 
rate. If (3 < 1, this process may be repeated to obtain 
the desired accuracy. For example, employing the second 
order form for the Newtonian potential 

?=-fSHO i +»»+4 (87) 



one can build the third order relation for the curvature 
$ 3 + 3w 



C 5 + 3w 



1 
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(88) 



In the opposite limit that the anisotropic stress con- 
tribution is large compared with the other perturbations 
(Sn > (, 5n > S), the integral solutions for the Newto- 
nian metric reduce to 



a J VP 
* = -$ + Sn ■ 



(89) 



Notice that the integrals remain finite as long as Sn di- 
verges at zero no faster than a ~ 5 / 2 ^ 3w / 2 anc l the prefac- 
tor yfp/a, is simply the decaying mode of the Newtonian 
curvature (see Eq. |59|). 

Applications. — Collisionless radiation provides an ap- 
plication for these results. The anisotropic stress of radi- 
ation is related to the quadrupole moment as defined in 
equation ([33|) by 



12^ 

n 7 = — e 2 



(90) 



on scales much smaller than the curvature radius. Equa- 
tion ( p3| ) also implicitly gives the equation of motion of 
the £th multipole as 



e, = k 
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(91) 



This is an infinite set of coupled differential equations 
representing the fact that radiative stress depends on in- 
ternal degrees of freedom of the distribution. However, 
|f42f introduced an approximation based on the solution 
of these equations in the absence of sources 



2£+l 

which allows us to express 



= Cje(krj) , 
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(92) 



(93) 



Applying this closure relation to the quadrupole, not- 
ing (a /a) — Iji] in the radiation-dominated era, and 
rewriting the dipolc in covariant form &i = v 1 — Ht 
shows that the anisotropic stress of free radiation in 
the radiation-dominated era obeys equation ( |si| ) with 

a = P-r/p = Py/ Pi-ad- 

The photons actually do not behave in this manner 
before last scattering since their coupling to the baryons 
destroys any quadrupole moment in the distribution (see 
j jVIIIE ). However, the same analysis applies to the 
massless neutrinos, whose anisotropic stress can be ap- 
proximately parameterized by a = Pv/p T &d- Equation 
( |S8| ) shows that the change in $ is enhanced by a fac- 
tor 1 + (2/15)(p I/ /p ra d). With equation (|56|), we find the 
effective temperature of the CMB to be 



e + * = = ~C(o) U-t^— 

2 3 \ 15 p rad 



and the ISW combination to be 
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(94) 



(95) 



to first order. These results are equivalent to but more 
physically transparent than those of pOj , In addition, 
equation (|8q ) introduces second and third order correc- 
tions that are important for models with additional neu- 
trino species or higher neutrino temperatures T v 



Pi 



0.681- 



N„ (1 M\T„ 



(96) 



as shown in Fig. 13 

Equations (|94|) and ( p5| ) imply that the neutrinos have 
significant but opposite effects on the effective temper- 
ature and ISW terms in the CMB anisotropy equa- 
tion d53|). The degeneracy is broken as the fluctuation 
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enters the horizon, as we shall see in § VIII D , leading to 
detectable effects in the acoustic peaks ]43(|. 

Finally, seed defect sources generally have large 
anisotropic stress contributions outside the horizon. 
Causality dictates that p s IT s behaves as white noise above 
the horizon which together with the so-called scaling 
ansatz leads to an r\~ x l 2 temporal behavior. This im- 
plies that Su oc ?7~ 5 / 2 oc a ~ 5 ( 1+3t1 ')/ 4 . Since this diverges 
slower than a -5/ 2 - 3lu / 2 j it does not imply a divergence 
of any observable |H| and contributes mainly to the de- 
caying mode of the curvature jjj. Note that these con- 
siderations assume vanishing spatial curvature. 



B. Entropic Stress 

We next consider the case in which entropic stress dom- 
inates the other stress components Se Ss,Su- In 
this limit Se does not depend on £ and equations ( |5l| ) 
and (|52l ) are more than simply a formal solution: they 
are the integral solutions for an arbitrary source in the 
absence of smooth components. For instance, if the en- 
tropic stress has a power-law behavior Se oc a n then 
the comoving curvature will have the same behavior 
C = — / dlnaSs oc a" assuming £(0) = 0. 

It is important to note that once Se turns off, £ will 
remain constant. Thus entropic stresses that act for only 
a fixed amount of time generate curvature fluctuations 
out of isocurvature initial conditions that then behave in 
the same manner as initial curvature fluctuations. 

A natural way of establishing an entropic stress is to 
have two components whose sound speeds differ at some 
point in the evolution. It is useful to introduce the "en- 
tropy" 



(5/92 



(97) 



P2 +P2 Pl+ Pi 

since its equation of motion is simply 

& = -k(v 2 - «i) (98) 

by virtue of equation ( jl4| ) , assuming no direct energy ex- 
change between species. The equivalence principle guar- 
antees that under purely gravitational evolution the ve- 
locity differences vanish so that a constant a is a good 
approximation outside the horizon. Tight coupling be- 
tween the components also implies a = 0. 

Let us assume that the density perturbations are ac- 
companied by sonic stresses in the rest frame of each 
component, 



A/o j = Spj - pj(vj -v)/k, 
Apj = Spj -pj{vj - v)/k, 



(99) 



with 



sumption vj 
sonic stress 



Apj/ Apj. Under the constant entropy as- 
v and with the definition of the combined 
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the entropic stress becomes 



St 



(9) 



(100) 
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Consequently, entropic stresses are generated when the 
sound speeds of the two components differ. Furthermore 
they can be much larger than the sonic stresses if Sp/ p <C 
a as is the case for isocurvature models. 

Full Solutions. — If the sonic stresses are also adi- 
abatic (c 2 = p'j/p'j), then Se = Sr and Sv/o = 
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(fJJp'y/3. Equation (|5l|) then yields 
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The Newtonian potential follows from equation (|5 
For example, in the constant w\ and W2 case 

Sp _ (l + wi)(l + w 2 )Aw 
a [1 + w\ + (1 + w 2 )y] 2y ' 
C = 1 (l + w 2 )y 
a 3 1 + wx + (1 + w 2 )y ' 



(103) 
(104) 



where recall Aw = w 2 — wi and y = p 2 / p\ and we have 
assumed isocurvature initial conditions (C(0) = 0). Once 
p 2 dominates the energy density, the entropic stress has 
been converted into a constant comoving curvature £ = 
cr/3. This solution may then be substituted into equation 
(|52|) to obtain the integral solution 
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where r = (1 + ioi)/(l + w 2 ) and n = — (5 + 3wi)/6Aw. 
The Newtonian curvature has the limits 
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The relation between the Newtonian and comoving cur- 
vature in the y 1 limit is exactly the same as the 



stress- free clustered case of §VI A 



Another interesting case is when the two components 
are initially both radiation iui(0) = w 2 (0) = 1/3 but the 
second component w 2 becomes non-relativistic. Equa- 
tion (102) implies that by the time p 2 3> pi, a curvature 
fluctuation 



a pi 
3 p 



(107) 
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is generated from isocurvature conditions. 

The more general solution (101) covers multicompo- 
nent models with non-adiabatic sonic stresses in one or 
both of the components and for example may be applied 
to scalar field models. 

Applications. — Entropy perturbations between the 
matter and the radiation are the basis of the prototypical 
isocurvature models, i.e. the baryon isocurvature model 
(PIB) and the axion isocurvature model. In this case, 
the integral in equation ( |105| ) can be explicitly solved to 
obtain [il 



$ _ 16 + 8y - 2y< 
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where y oc a. Thus, the curvature grows as a in the 
radiation-dominated era only to freeze out at the ampli- 
tude reached at matter-radiation equality. 

Motivated by our study, one of us constructed an 
isocurvature decaying dark matter (iDDM) model that 
utilizes the mechanism described above where density 
fluctuations in two radiation species are initially balanced 
p2| . The resulting constant superhorizon curvature fluc- 
tuation £ is a property generally associated with adia- 
batic models. Another interesting consequence is that 
the ratio of large-scale CMB temperature anisotropics to 
the Newtonian potential is neither 1/3 nor 2. This is 
because the photons possess initial perturbations. These 
can be chosen to make the ratio greater than or less that 
2; in particular it can be arranged to be close to the 
adiabatic 1/3 relation. 



C. Sonic Stress 



Sonic stresses provide the final pure case. Here 
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P + P 
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recall that adiabatic stresses are a special case of a sonic 
stress where c 2 c — c 2 s = p'/p'. Because the stress fluc- 
tuation is related to the comoving density fluctuation by 



equation (109) and in turn to the Newtonian curvature 
via the hybrid Poisson equation (p4|) , the evolution equa- 



tion (50) for the Newtonian curvature becomes a homo- 



geneous second order differential equation. 
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dp, 
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where the sound horizon is defined as 



dijee ■ 



(111) 



Here, we have assumed that the wavelength is much 
smaller than the curvature scale, but all results are ap- 
plicable to the general case with the replacement k — > 
Vk 2 - 3K. 

For ks' 3> 1, we can approximate this as an oscillator 
equation with an effective mass 



fcV 2 m eff $ = , 



with 
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Under the WKB assumption that the effective mass is 
slowly-varying compared with the frequency of oscilla- 
tion, the fundamental solutions to this equation are 



cc 
cc 
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cos(fcs) , 
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sin(fcs) 



(114) 



This equation says that once the fluctuation passes inside 
the sound horizon, the Newtonian potential oscillates and 
decays reflecting analogous behavior in the density per- 
turbation through the Poisson equation. 

Full Solutions. - - Equation (114) suggests that a 
change of variables to 
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cc 
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should simplify the equations of motion. Indeed, equa- 



tion (11C) becomes 



where 
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If F is constant, then equation (116) is a variant of 
the Bessel equation, and the solutions Q will approach 
sin(/cs) (up to a phase) at large ks. Furthermore, when 
ks ^ the term with F may simply be neglected, 

and again the solutions will be sin(/cs) with arbitrary 
phase. Because the latter solutions will match trivially 
onto the former, we reach the conclusion that if, for 
a given mode, F is constant until s 3> \F\ x / 2 /k, then 
the appropriate Bessel function solution will hold for all 
times, regardless of how F varies once ks S> l-Fl 1 / 2 . If 
cc and p vary on the Hubble time scale, then F is or- 
der unity, and the solution describes modes that are well 
inside the horizon before F begins to vary. 
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The simplest way to arrange F to be constant is to have 
a constant equation of state w\ and constant sound speed 
cc- For constant 101, s is only defined for w\ > —1/3, 
and we therefore set p — p cr . Then 



F = 6 



1 + Wi 



(118) 



independent of the sound speed. The growing and de- 
caying solutions for <I> arc then 
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(119) 



where v = (5 + 3u>i)/(2 + Qw\). Again, once a given 
mode is well inside the sound horizon, the condition that 
wi and cc be constant can be relaxed. 

Applications. — Adiabatic stress perturbations in a 
baryon-photon universe provide the prototypical exam- 
ple. In this case, the tight coupling between the photons 
and baryons through Compton scattering prevents the 
generation of entropy through the motion of matter (see 



VIII C). The growing mode of equation (Ul9) becomes 
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with k y = (a/a) y= i and 
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We have chosen the normalization to match the stress- 
free solution <E> = 2^/3 initially. In this case, the acoustic 
oscillations of the photons are directly related to behavior 
of the potential via the Poisson equation 



v 



1 + 3y/4 



(l + 3y/4)- 1/4 C(0)cos(fcs) (fcs>l). (122) 



Since the acoustic oscillations are responsible for the 
acoustic peaks in the CMB from equation (]5^) , this deter- 
mines the morphology of the acoustic peaks in this simple 
adiabatic photon-baryon universe. The important result 
is that the acoustic oscillations follow a cos(fcs) pattern 
in phase with an amplitude that is enhanced by the decay 
of the initial curvature perturbation. These results are in 
fact generic to adiabatic models due to similar evolution 
in the driving potentials [ OJ , 

A second example is provided by the QCDM and 
GDM models where the equation of state can go neg- 
ative (w < 0) while the comoving sound speed remains 
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FIG. 14. Sonic stresses in a scalar-field dominated universe. 
Sonic stresses cause scalar-field perturbations to stop growing 
and hence lead to a decay in the Newtonian curvature in a 
scalar-field dominated universe. Here we compare a QCDM 
model with wq — and no CDM with the analytic prediction 
of equation (119). The discrepancy at early times is due to 
radiation contributions in the QCDM model and those at late 
times from the baryons. 



real (c c > 0). In QCDM, c c = 1 by virtue of the scalar- 
field equations of motion. In GDM, it is allowed to have 
any positive value. The solutions are strictly valid only 
after the exotic component dominates the expansion rate 
and fluctuations. To test this solution, we show an ex- 
treme example where the dark matter is all in a "Q" 
component (no CDM) with wq = (see Fig. |l4|). The 
small departures at early and late times are due to the 
radiation and baryonic contributions. 



VIII. MIXED STRESSES 

Although the purely anisotropic, entropic and sonic 
stress cases of the last section illustrate many aspects of 
stress phenomenology, when all three are present they 
can interact and create a diverse range of behavior. In 
this section, we study a few typical cases: sonic stress 
generation from entropic stress, entropic stress genera- 
tion from sonic stress, and anisotropic stress dissipation 
of sonic stress. The first process is responsible for gen- 
erating acoustic phenomena in all isocurvature models. 
The second process is responsible for the heat conduction 
in fluids and the generation of smooth components. The 
last process is responsible for viscous dissipation of acous- 
tic phenomena and also provides an alternate means of 
generating a smooth component. We briefly discuss seed 
(defect) stresses, which carry not only all three types of 
scalar stress but also vector and tensor stresses as well. 
Finally, we consider the effect of tensor anisotropic stress 
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on gravitational waves and CMB anisotropics in passive 
models. 

A. Formal Solution 

If the equations of motion for a sonic stress pertur- 
bation are known to be solved by $1 and $2, the full 
solution can be written in integral form as, 



ra da $ 1 (d)$ 2 (a) - $i(a)$ 2 (a) 

1 P'cr + P' _ P 

2 Per P 



(123) 
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where A and B are arbitrary constants. This solution 
is merely formal unless the remaining stress sources can 
be specified independently of We now consider some 
special cases. 



B. Entropic/Sonic Stress 

A model that begins with isocurvature initial condi- 
tions will generate adiabatic density fluctuations through 
the action of entropic stresses, as discussed in § VII B . 



These density perturbations carry with them sonic per- 
turbations through equation (|3^) that stop the further 
growth of density perturbations throu gh the same pres- 
sure support mechanism discussed in § VII C . Therefore, 
entropic growth of curvature fluctuations generally ceases 
once the fluctuations cross the sound horizon of the dom- 
inant species. 

Full Solutions. — The solutions for $ in equati on ( 11£) 



along with the entropic stress from equation (101) in 
equation (123) allow us to construct the full solution in 
the presence of a constant entropy and p = p CI 
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Recall that the combined and component sound speeds 
were defined in equations (10C) and ( |99| ) respectively. 
This solution is valid for k/k y ^> 1 and k ^> IK] 1 / 2 (see 

sfvnc]). 

Applications. — The case of ui\ — 1/3 and w 2 = 
is of special interest because it corresponds to baryon- 
photon entropy perturbations, as in the PIB model before 
last scattering. The result is that on small scales where 
ks S> 1 and k/k y >• 1, the curvature behaves as |3|] 
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which is again direct ly re lated to the temperature oscil- 
lations by equation (122). The first term in brackets is 
due to density perturbations in the baryons remaining 
from the constant entropy condition cr = 6b — §<5 7 ~ St,. 
The second term represents decaying acoustic oscillations 
from the adiabatic pressure. The extra factor of k y /k re- 
flects the fact that the curvature grows as a until sound 
horizon crossing at an oc (k y /k) Q. 

An interesting result is that the acoustic oscillation fol- 
lows a sin(fcs) relation, implying the opposite phase in the 
acoustic peaks compared with the cos(fcs) adiabatic case 
]l0t . This result is rather generic to isocurvature models 
again due to the similar behavior of the driving poten- 
tials. For example, axionic isocurvature models where 
the entropy is between the radiation and the CDM also 
follows this pattern. 

For this reason, isocurvature seed pressure also tends to 
generate this type of acoustic pattern |ll| . These stresses 
are found in topological defect models; indeed the dom- 
inant scalar mode of strings, monopoles, textures, and 
non-topological textures do behave in this manner p3[ . 
However, only in the latter two are the other modes suffi- 
ciently small at the first few peaks to yield clean acoustic 
peaks even for the scalar perturbations alone. Defect 
models generically have vector and tensor stresses that 
generate comparable levels of anisotropy and further ob- 
scure acoustic phenomena. 

It is possible to construct isocurvature models with 
an "adiabatic" pattern of acoustic peaks. The simplest 
way to arrange this is to create constant comoving curva- 
ture perturbations £ through entropic stresses that turn 
off well before the perturbation crosses the horizon. A 
concrete example of this kind of mechanism is given in 
I . As is clear from § VII B and originally pointed out by 
11 1, for an isocurvature model to generate scale-invariant 
curvature perturbations requires that the entropic stress 
have superhorizon scale correlations which cannot be gen- 
erated causally without an inflationary epoch. 



C. Sonic/Entropic Stress 

Likewise, an adiabatic or sonic fluctuation will not gen- 
erally remain so as it crosses the horizon. Inside the hori- 
zon, the fact that components with different equations 
of state have different pressure responses to gravitational 
compression will cause the species to move independently. 
The generation of entropic stresses is in fact a primary 
mechanism for creating the smooth components of mat- 
ter discussed in 



VI B 



Full Solutions. — The case of two components with 
constant ui\ and wi again provides an instructive ex- 
ample. The equation of motion for the entropy (p8) is 
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constructed out of combining the continuity equations 
(U) of the two species and has the formal solution 



(127) 



a = er(0) — k / drj(v2 ~ V\) 



Recall th e en tropy is related to the entropic stress by 
equation (101). Note that entropy leads to entropic stress 
only if the two components differ in their equation of 
states (Aw ^ 0). 

Since we are interested in the generation of entropy, 
we will assume that the initial entropy perturbation <r(0) 
vanishes. Two interesting cases are when the entropy 
generation timescale 7?r is small compared to the expan- 
sion time and when it is comparable to the expansion 
timescale. 

In the former case, the entropy will lead to dissipative 
behavior if 



a w —krjrv , 



(128) 



where rjr is some characteristic timescale for entropy gen- 
eration. While this form is not gauge invariant even 
though a is, the ambiguity vanishes inside the horizon. 
Hence, equation ( |128| ) is a good approximation in the 
desired case rjr <§C r\. Substituting equation fll28p into 
(101) and then (Jl^) in Newtonian gauge and assuming 
a solution of the form v oc exp(i J ujdrj), we obtain the 
dispersion relation 
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where recall cc was defined in equation QlOOj) . Note that 
we have assumed rjr/rj <C 1 in order to replace rest-frame 
sound speeds with comoving sound speeds. 

This describes an exponential damping of sound waves 
as the wavelength passes the "diffusion" scale k — 

In the opposite regime, where the entropy is generated 
on the horizon scale, the sonic nature of the total system 
breaks down before acoustic oscillations even start. In 
this case, components can decouple completely and form 
subsystems where sonic, anisotropic, or smooth effects 
can occur separately. Here it is simpler to describe the 
behavior of individual subsystems through variables that 
are comoving with respect to the individual species J. 
The sound speed cj and the contribution to the Newto- 
nian potential are useful (see Eq. p9]) 
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which implies <f> = ^2 j If the density fluctuations in 
the other species are still below their own sound horizon, 
then the oscillating components can become effectively 
smooth in comparison. From that time forward, we can 
treat the system as having a smooth component, and the 




FIG. 15. Creation of a smooth radiation background. Nu- 
merical results in a ACDM model are shown here. Notice 
that upon crossing the sound horizon, the photon contribution 
$ 7 to the total curvature $ damps and oscillates, while the 
neutrino contributions damp much more rapidly due to colli- 
sionless damping from its anisotropic stress. The oscillating 
photon perturbations yield little time-averaged effect, and the 
CDM evolves under equation ([n]) leading to the well-known 
logarithmic tail in the CDM transfer function |ft5j. 



curvature fluctuation contributed by the other species is 
governed by equation j69|). 

Applications. — The case of a short entropic timescale 
Vt *C t) is realized in the photon-baryon fluid before re- 
combination and is relevant for considering the damp- 
ing of acoustic phenomena in the CMB for all models of 
structure formation Jl6],fl7) . Here the entropic time scale 
is derived from the baryon Euler equation 
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where R = Spb/^p-y- Inserting this into equation (127) 
we obtain rjr — t~ 1 R under the rapid oscillation assump- 
tion. The dispersion relation then becomes (Iffl 



w=±fecs+ 4 fcV " 1 (iw 
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These entropic pressure techniques provide a simpler and 
more transparent derivation of this well-known result 
than exists in the literature. Note that heat-conduction 
damping is suppressed by R in the photon-dominated 
epoch, and we shall see that in that case viscous damp- 
ing from the anisotropic stress is more important. 

The opposite limit of an entropic timescale on order 
the expansion time rj-p ~ rj is applicable to all models 
with CDM. The CDM never participates in the acoustic 
oscillations of the baryon-photon system even during the 
radiation-dominated era. In Fig. we show an example 
of a perturbation deep in the radiation-dominated era in 
a CDM model with the usual neutrino content. Before 



24 



horizon crossing kr\ <C 1, the perturbations are adiabatic 
and the total potential $ is constant with contributions 
from the photons <i> 7 and the neutrinos <&„. The sound 
horizon for the radiation is s = ry/VS and once kr\ ^> \/3, 
the potential contribution of the photons starts to decay 
as in equation (12C ) . The neutrino contributions decay 
even faster due to their anisotropic stress, as we will dis- 
cuss in the next section. The CDM contribution $cdm 
then turns over and behaves as 4>cdm oc ln(Ca) / a where 
C is some constant This is exactly the behavior 



predicted by equation (|69D assuming the radiation com- 
ponent is smooth. In actuality, the photon component 
is not smooth compared with the CDM component but 
oscillates sufficiently rapidly that its time-average is neg- 
ligible. 

These considerations also apply to models in which an 
additional component with iogdm < comes to domi- 
nate at late times. Quintessence a nd GD M models (see 
Fig. [l6|) are examples thereof. In § VII C , we considered 
the case where the GDM completely dominated the ex- 
pansion. Here we consider the effect of adding a compo- 
nent of CDM. 

The critical parameter here is the comoving sound 
speed of the GDM cqdm- Slowly-roling scalar fields 
found in the QCDM subcategory have cqdm = 1, but 
components with lower sound speeds are possible in prin- 
ciple. The sound speed tells us how long after horizon 
crossing the GDM component stabilizes due to pressure 
support. If the perturbation is already within the sound 
horizon by the time the GDM comes to dominate the ex- 
pansion rate (y = 1), then the total potential will behave 
under equation (|6^) as if the potential is smooth regard- 
less of the exact value of the sound speed. As the sound 
speed is lowered such that crossing occurs near (y = 1), 
we see effects from the finite sound speed. In the limit 
where the perturbation remains above the sound horizon 
until the present, the solution returns to the clustered 
case of equation (M). 



D. Sonic/ Anisotropic 



We have seen in § VII A how anisotropic stress can gen- 
erate Newtonian curvature perturbations but, like adia- 
batic stress, it can also destroy them. In this context, 
anisotropic stress represents the "frictional force" set up 
in response to the non-uniform bulk flow of matter and 
shear in the metric. As we shall see, this sort of behavior 
is not confined to fluids. It represents another mechanism 
for generating a smooth component. 

Full Solutions. — The phenomenological parameteriza- 
tion of anisotropic stress in equation (|8^) yields two in- 
teresting limits. The short timescale limit (kr/u <C 1, 
Vu/v 1) leads to viscous or collisional damping. Here, 
the anisotropic stress is algebraically related to the ve- 
locity in shear-free frames 




FIG. 16. Creation of a smooth quintessence or GDM com- 
ponent with icqdm = —1/3. As shown in Fig. |l4|, sonic 
stresses in the GDM component prevent perturbation growth 
inside the sound horizon CGoyikn ~ 1. If this occurs well be- 
fore GDM domination at y = 1, then the total Newtonian 
curvature will evolve as if the component were always smooth 
(i.e. under equation |7l| ). If this occurs well after, the New- 
tonian curvature will behave as if all components were fully 
clustered (i.e. under equation [^3|). For a given scale, this 
depends on the comoving sound speed cgdm; for quintessence 
cgdm = 1. 



n = 4ar)n(kv - H T ) . 
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Under the same assumptions used to derive the disper- 
sion relation for heat conduction ( |12S| ), we find the vis- 
cous dispersion relation 
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implying dissipation at a characteristic scale k ~ 

If the dissipation timescale is comparable to the expan- 
sion timescale, the damping occurs at horizon crossing 
but is more gradual. The formal solution in this limit 
is given in equation (|82|) and approximates the effects of 
collisionless damping. 

Applications. — Collisional damping occurs in the 
photon-baryon fluid before last scattering and is the pri- 
mary dissipation mechanism for acoustic oscillations the 
CMB jl6| . Equation ( |133| ) then describes the anisotropic 
stress of the photons II 7 with a = 2/5 and ryn 
Repeating the calculation leading to equation (J13S 
dispersion relation for the oscillations becomes ]46|| 



uj = ±kcc + i 
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In comparison to the heat conduction dissipation of equa- 
tion ( 132 ) , viscous dissipation is more effective if R < 1 , 
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which is the case for the baryon content implied by big 
bang nucleosynthesis (e.g. 

Collisionless damping occurs in free radiation. Free ra- 
diation behaves in this "frictional" manner because gra- 
dients in the potential flow are dissipated as radiation 
streams from one part of the flow to the next. In terms 
of the multipole moments, power in the dipole gets trans- 
ferred to the quadrupole and so on through the hierarchy 
equations ([H]). As such, the anisotropic stress acts as the 
gateway for anisotropy generation in the CMB. It also 
leads to a damping of density and velocity perturbations 
inside the horizon for all species of free radiation and 
is the reason that in Fig. [IB] the neutrino contributions 
damp more rapidly than the photon contributions. More 
generally, it is responsible for smoothing out components 
when entropic stresses cannot be generated, e.g. when 

c|-C? = 0. _ 

These dissipational terms may also be important in 
stabilizing other forms of matter. One way to generate 
a smooth density component is to introduce an effective 
viscous rather than sonic stress [|l7| ; a mechanism of this 
type is thought to be involved in stabilizing the strCDM 
model [^6|. Turok ]5(J] suggested an extreme example 
of this sort, in which the comoving sound speed of the 
seeds is imaginary but the anisotropic stress is perfectly 
balanced to counter the otherwise exponential growth of 
density fluctuations. 



E. Sonic/Entropic/ Anisotropic Seed Stress 

The seed stresses of topological defect models provide 
an example where all types of stresses coexist. While the 
consequences of a given seed stress for structure forma- 
tion are straightforward to work out, the behavior of the 
seed stress itself is more difficult and requires simulations 
to work out in detail. 

Full Solutions. — The two-point statistics of compli- 
cated defect models can be accurately modeled as the 
incoherent (quadrature) sum of a relatively small num- 
ber (~ 10-20) of simple seed stress histories pqjl3fl . Each 
individual source may then be determined by the tech- 
niques above (see Eq. |Q and Ea| ). 

Applications. — The simplest defect models typically 
have two other properties that have important phe- 
nomenological consequences. Defect models are causal in 
the classical sense. The stress perturbations <5p s and p B Tl s 
must fall off at least as white noise (k°) outside the hori- 
zon and the initial curvature must vanish (£(0) = 0). The 
traditional string and texture models also obey a scal- 
ing relation that states that the stress histories depend 
on wavelength only through the combination kr\. This 
ensures self-similarity of the structure at horizon cross- 
ing and leads to nearly scale- invariant CMB anisotropies. 
Because the simplest versions of these models run into 
difficulties when CMB anisotropies are compared with 



large-scale structure, phenomenological models that do 
not obey the scaling relation have recently received some 
attention [E5l. 



F. Vector Stress 

We next consider the effect of vector anisotropic 
stresses. Recall that in the absence of vector stress, 
the vector perturbation decays. In order to generate an 
observable effect, vector perturbations must be continu- 
ously generated by vector anisotropic stresses. This im- 
plies that vector modes are unique to the active stress 
models, of which defects are an example. 

Full Solutions. — The solutions in the presence of stress 
sources can be constructed via Green's function tech- 
niques from the stress- free solutions of (|7^). For the 
vector modes, the solution becomes 
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where C =const. represents the decaying mode. The re- 
maining metric perturbation kB^ 1 ^ — H^ 1 ^ is related 



algebraically to equation (136) through equation (pT|). 
Applications. — The discussion of scalar seed stress in 



defect models in § VIII E also apply to vector stresses with 
a few additional considerations. Defect models generally 
have comparable scalar, vector, and tensor anisotropic 
stress sources above the horizon f51|| . Since CMB 
anisotropies are primarily generated at horizon crossing, 
these sources tend to yield comparable anisotropy contri- 
butions for modes that cross after last scattering. Vector 
modes that cross before last scattering do not contribute 
due to suppression of tensor anisotropies in the CMB 
from scattering. Thus, vector modes can obscure the 
first few scalar acoustic features in defect models. On 
the other hand, vector modes have a special signature in 
the CMB polarization that may assist in their isolation 



G. Tensor Stresses 

Tensor anisotropic stresses provide sources and sinks 
for gravity waves. It is well known that quadrupolar 
stresses generate gravity waves. Furthermore, a passing 
gravity wave will also impart some energy to the radi- 
ation backgrounds via differential gravitational redshifts 
and thereby decay. 

Full Solutions. — An integral solution may be con- 
structed out of the homogeneous solutions of § VI D as 
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H, 



(±2; 



{r,) = c 1 H 1 {< n ) + c 2 H 2 {n) 



(137) 



^ ^ H 1 (f l )H 2 (f 1 )-H 1 (r 1 )H 2 (f l ) o 2 (±2) 
V^^-^W)^^ 8 ^ PH ' 



where Ci and C2 are constants associated with the initial 
conditions. 

Applications. — We show the damping effect of 
anisotropic stress in the radiation backgrounds in Fig. [Tt]. 
It is generally not included in standard Boltzmann codes 
that solve CMB anisotropics |54j . The reason is that it is 
negligible at large angles since the corresponding modes 
entered the horizon well into matter domination when the 
anisotropic stresses are negligible. Since this feedback 
effect is typically small, equation (137) may in princi- 



ple be used to iterate to a solution from the free-gravity 
wave case of equation (|r8|). However, the evolution of the 
anisotropic stress has a more important effect. Just as 
scalar anisotropic stresses in the photons are destroyed 



by scattering before last scattering (see §VIIID), tensor 
anisotropic stresses are destroyed as e~ T . This cuts off 
the tensor contributions to the anisotropics before last 
scattering as indicated in equation (|7^) and shown in 
Fig. [l7|. The feedback effect still exists, but the level of 
the anisotropy itself makes it too low to be observable for 
reasonable tensor to scalar ratios. 



The Green's function solution of equation (137) is 
more directly applicable to the seed stress case where 
pH ±2 — p s Il s . As in the vector stress case, the phe- 
nomenological result is that defect models tend to have 
significant tensor contributions above the angle the hori- 
zon subtends at last scattering. As in the passive models, 
their contribution is cut off below this scale due to scat- 
tering. 



IX. A DESIGNER APPLICATION 

With this general study in hand, we are now in a po- 
sition to discuss prospects for reverse engineering the 
model for structure formation. Obviously, the specific 
route the inversion takes will depend on the results of 
ongoing experiments. Currently, the data favor a model 
with phenomenology like the ACDM model, e.g. the 
shape of the large-scale structure power spectrum |p5,p6|, 
relative high to low redshift cluster abundances |p7,p8|, 
supernova luminosity distances to redshift of z ~ 0.5 
p9|-|6l]|, and degree scale CMB anisotropies (e.g. ]62| , |63] ] ) . 
If agreement between the ACDM model and future preci- 
sion tests is good, can we say purely on phenomenological 
grounds that we have proven the existence of a cosmo- 
logical constant? If the model varies from the data, can 
we use the methods developed here to modify the stress 
history, restore agreement, and in the process learn new 
information about the dark sector? 

Let us address the first question. We know that the 
evolution of structure is completely defined by the stress 




FIG. 17. The effect of tensor anisotropic stress on tensor 
anisotropies. Gravity waves generate tensor anisotropic stress 
in radiation that absorbs energy and damps the gravity wave. 
For free radiation like the neutrino background radiation, this 
reduces the anisotropies generated on small scales. For the 
CMB, the elimination of tensor anisotropic stresses by scat- 
tering cuts off small-scale contributions as well. 

history of the matter. Since the equations of state of 
the ordinary matter are known, the remaining element 
is the dark sector. To test the uniqueness of the ACDM 
model, we should look for alternate means of reproducing 
its stress history. The background stress history of the 
dark sector in a A model is given by 



WACDM 



a 3 + ^cdm/^a ' 



(138) 



and its stress perturbations vanish. This suggests that a 
generalized dark matter (GDM) component of the type 
introduced by (lj] and parameters 



WGDM — WACDM, 



C GDM — I 



(139) 



should reproduce the phenomenology of the ACDM 
model. Recall that cqdm is the s ound sp eed in the frame 
comoving with the GDM (see § |VIII C| ). To the extent 
that the comoving and GDM-comoving frames coincide, 
this form of dark matter exactly reproduces any mixture 
of A and CDM. We show an example in Figure [H| Note 
that all classical cosmological and linear theory tests will 
return the same answer for the two models despite the 
fact that the GDM model is a single dark matter com- 
ponent model in a critical density universe! Non-linear 
effects are also identical if the same stress history is main- 
tained throughout. 

There are two lessons here. The first is that on purely 
phenomenological grounds we can do no more than mea- 
sure the global properties of the dark sector. A multi- 
component model and a single component model with 
the same stress history are formally identical. 

The second is that reproducing the phenomenology of a 
ACDM model is rather simple: it requires an equation of 
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FIG. 18. Mimicking ACDM. Models with the same stress 
histories in the dark sector have the same observable conse- 
quences regardless of differences in how that sector is com- 
posed. Here numerical solutions for the CMB anisotropy and 
large-scale structure power spectra in a ACDM universe and a 
single -com ponent critical GDM dominated universe of equa- 
tion (139) are shown to be identical. However, variations in 
cc, A$, and Arj away from this tuned stress history do have 
observable consequences. 

state that varies from w = at high redshift to w w —2/3 
today and a form of matter that is free of large-scale 
stress gradients in its comoving frame. Thus, there exists 
a wide class of such single component models that fit the 
current data as well as the ACDM model. 

On the oth er hand, variations on the conditions in 
equation ( 13£ ) have observable consequences for future 
measurements. Relaxing the stress-free perturbation 
condition by raising the sound speed reduces the small- 
scale power in the model and delays the formation of high 
redshift objects. The remaining freedom in the stress his- 
tory is associated with the equation of state w. We can 
quantify this by recalling that the gravitational potential 
depends only on the quantity (see equation |p9[ ) 

A$ oc & f * . (140) 
a Jo VP 

Large-scale structure constrains the value of this inte- 
gral at the present (a = 1). CMB anisotropics from the 
ISW effect are sensitive to variations in this function that 
occur on the order of the light-travel time across a wave- 
length (see §VD). Since CMB anisotropies potentially 
probe nearly three orders of magnitude from the current 
horizon, variations at a fraction of a percent of the cur- 
rent expansion time are potentially visible in the CMB. 

Similarly, distance measures such as the angular diam- 
eter distance to the last scattering surface (a ~ 10 -3 ) and 



the luminosity distance to high redshift objects (a 
probe the combination 



Arj(a) oc 



da 



1 



aV /2 ' 



2/3) 



(141) 



Current measures of A$(a) and A-q(a) are crude at 
best as they only proble their values at discrete epochs 
or averaged over long timescales. A sharp test of the 
A+CDM hypothesis that should be possible with future 
measures involves reconstructing the time evolution of 
the equation of state u> cr through measures of A<I>(a) 
and/or Arj(a). Any combination of a cosmological con- 
stant and CDM will obey 



3w cr (l + w cr ) 



(142) 



The physical implication of this relation is that the back- 
ground pressure is constant in time (p' = 0) so that the 
adiabatic sound speed vanishes (c^ = p'/p' = 0). If this 
condition is violated, we will have proven the existence 
of a new form of matter. If w' CI > 3ui cr (l + iv CI ), then it 
can be supported by adiabatic stresses f64j] : if not, then 
a form of matter with non-adiabatic stresses is required. 
Non-adiabatic stresses imply that the relation between 
the pressure and density perturbations does not follow 
that of their spatial averages. 

The simplicity of the requirements of equation ( |l39| ) 
raises the possibility of a unified description of the dark 
sector. A concrete but somewhat trivial example is a 
scalar field that rapidly oscillates around a non-zero po- 
tential minimum. The rapid oscillations average away all 
large-scale pressure effects save from the vacuum pres- 
sure of the potential minimum J65| ] . Unfortunately, the 
relationship between the constant and quadratic pieces of 
the scalar potential is left unexplained and so is no bet- 
ter than an explicit A+CDM model. Nonetheless, there 
may be more complicated examples, perhaps involving 
multiple fields, in which the mimic conditions ( |139| ) are 
approximately satisfied and that do unify the two behav- 
iors in a true sense. 

This discussion shows that a reverse-engineered model 
for structure formation will in general not be unique. On 
the other hand, the observables can be translated into 
constraints on the stress histories and phenomenological 
models of the dark sector. These in turn can assist in the 
search for compelling physical candidates to compose the 
dark sector. 



X. DISCUSSION 

Without any assumptions besides the validity of gen- 
eral relativity and nearly homogeneous and isotropic ini- 
tial conditions, the evolution of structure is completely 
determined by the stress history of matter. We have stud- 
ied the means by which stresses, both in the background 
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and the fluctuations, can alter the observable properties 
of the model. 

We have examined the effects of smooth, anisotropic, 
sonic, and entropic stresses in structure formation, in- 
cluding their interactions and ability to generate effec- 
tively smooth density components. We have illustrated 
these behaviors with analytic solutions for systems with 
multiple components of differing background equations 
of state, which can themselves be time dependent in sev- 
eral important cases. These solutions have applications 
to nearly all of the current models for structure formation 
and are substantially more general than those existing in 
the literature. 

Although this study is not exhaustive, we have made 
explicit all of the assumption required to arrive at specific 
models and their accompanying phenomenology. In the 
process, we have exposed the limitations of traditional 
categorization schemes like that in Fig. [j]. These distinc- 
tions can be blurred in cases where the usual assumptions 
do not apply. We summarize several notable cases here: 

Initial Conditions: 

Isocurvature initial conditions imply a growing comov- 
ing curvature outside the horizon on scales relevant to 
large-scale structure and degree-scale anisotropics. The 
comoving curvature grows outside the horizon only by 
the action of stress perturbations. Once stress pertur- 
bations are turned off, the curvature remains constant 
until horizon crossing or curvature domination. These 
considerations provide a means for mimicking the phe- 
nomenology of adiabatic models |^2| and are important 
for interpreting the implications of CMB acoustic peak 
phenomenology; however, they do not invalidate the con- 
clusions of pd| ]. 

The Newtonian curvature $ is simply proportional to 
the comoving curvature if the background equation of 
state is constant. The Newtonian curvature admits a de- 
caying mode whereas the comoving curvature does not. 
The decaying mode can be stimulated by anisotropic 
stress perturbations outside the horizon but has observ- 
able consequences only through the contribution remain- 
ing at horizon crossing . 

// the Newtonian potential ^ is constant from last scat- 
tering to the present, the observed temperature perturba- 
tion depends only on the equation of state and ^ . The 
assumption here is that the comoving temperature per- 
turbation is negligible and is only true if stress pertur- 
bations are also negligible compared with the comoving 
curvature for all time. The axion isocurvature model pro- 
vides a simple counterexample. On the other hand, no 
assumptions about the anisotropic stress are necessary 
even when $ = no longer holds. 

Isocurvature initial conditions predict an observed tem- 
perature perturbation of 2^ on scales larger than the 
horizon at last scattering. The assumption here is that 



the initial temperature perturbation in Newtonian gauge 
vanishes. This is not the case for models where the 
isocurvature conditions are established by balancing the 
photon density perturbations off another species of ra- 
diation. Furthermore, changes in the potential that are 
slowly varying compared with the light travel time across 
a perturbation do not affect the observed temperature. 

Clustering Properties: 

The smoothness of a component is gauge-dependent 
and hence has no physical meaning. The gauge depen- 
dence of a smooth component is not a problem per se 
as certain frames, e.g. the frame where the momentum 
of the component vanishs, are dynamically special. A 
smooth contribution to the density with ws ^ — 1 does 
violate covariant energy conservation in any coordinate 
system where the spatial curvature changes. Since the 
very presence of a smooth density component requires 
that the comoving curvature perturbation decays, there 
can be no identically smooth contributions in those coor- 
dinates except in the trivial zero curvature perturbation 
case. A component can be smooth relative to another 
species inside the horizon where the relativistic effects of 
curvature variation are negligible fl5| . 

The behavior of a smooth component depends only on 
its equation of state ws ■ Since all components except A 
and curvature are clustered outside the horizon, the man- 
ner in which a component becomes smooth is observable. 
The presence of sonic (supportive) and anisotropic (dis- 
sipative) stresses are two possibilities, but others exist in 
principle. 

A "clustered" model like standard CD M has no dynam- 
ically important smooth stresses or density contributions. 
The radiation backgrounds and baryons are effectively 
smooth well inside the horizon prior to recombination. 
Smooth contributions are generic to structure formation 
models. 

Smooth component behavior implies small density per- 
turbations dps < 5 pc- A component can be effectively 
smooth even while possessing large density fluctuations. 
The crucial assumption is that their time-average den- 
sity over the dynamical time of the clustered species 
is smooth. Density perturbations in a component that 
vary rapidly with the expansion time generally lead to 
no effect on the growth of structure. The radiation back- 
grounds mentioned above provide a familiar example. 

The missing mass (clustered dark matter) and miss- 
ing energy (smooth ws < —1/3 dark matter) are sep- 
arate problems. The stress history of the dark sector 
completely defines its properties for classical cosmology 
and structure formation. Any combination of compo- 
nents that produces the same stress history will produce 
the same phenomenology. As an example, we have con- 
structed a toy model that that exactly reproduces the 
ACDM phenomenology but employs a single component 
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of dark matter in a critical density universe. Variations in 
the stress history produce models that satisfy the current 
constraints equally well but are potentially distinguish- 
able from ACDM. 

Perturbation Type: 

Scale-invariant gravity waves preferentially enhance 
the low- order multipoles of CMB anisotropics. Enhance- 
ment only occurs if the gravity wave amplitude decays 
close to horizon crossing and is eliminated as the equa- 
tion of state of the background drops. 

There is a sharp distinction between active and passive 
models for structure formation. Models that have been 
labeled "passive" in the literature are those in which the 
stress perturbations are simply related to density, veloc- 
ity, and metric fluctuations by equations of state, sound 
speeds, viscosity parameters, etc. Models in which there 
is a component whose stresses have no fixed relation to 
density and metric perturbations have been labeled "ac- 
tive" . The issue is the number of internal degrees of 
freedom that act to specify the stresses. For example, 
scaling defect stress histories are typically approximated 
by tens of parameters and those of particle dark matter 
by three or fewer. In principle, there is a spectrum of pos- 
sibilities between the two and a corresponding spectrum 
of phenomenological consequences. 

Our study is useful even if the current evidence sup- 
porting ACDM-type phenomenology holds up. Even the 
ACDM model itself does not have an entirely trivial stress 
history as its dark sector includes the neutrino back- 
ground radiation. The observability of the neutrino stress 
history has been addressed numerically in |J43| ; we have 
examined its physical origin here. Furthermore, the dark 
sector could contain exotic features that produce more or 
less ordinary phenomenology, and one needs to construct 
sharp tests against alternatives. For example, a combi- 
nation of cold dark matter and a cosmological constant 
must obey w' CT = 3u> cr (l + w cr ). This relation also acts 
as the dividing line between models with exotic and or- 
dinary matter. For exotic matter, the stress and density 
perturbations obey a relation that opposes that between 
the background stress and density; scalar fields are one 
example of exotic matter. 

In summary, a purely phenomenological reverse- 
engineering of the model for structure formation will re- 
quire the reconstruction of the time-averaged stress his- 
tory of the dark sector. This inversion is generally not 
unique. Nonetheless, if the observed phenomenology re- 
mains close to that of our simplest models, our study 
of stress phenomenology should provide the means for 
constructing viable models. If the observations require 
more radical departures, our study should be useful in 
identifying the assumptions that are incorrect and assist 
in the search for the correct phenomenological model for 



structure formation. 
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